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NATIONAL  BUREAU  OF  STANDARDS 


The  National  Bureau  of  Standards1  was  established  by  an  act  of  Congress  March 
3,  1901.  Today,  in  addition  to  serving  as  the  Nation's  central  measurement 
laboratory,  the  Bureau  is  a  principal  focal  point  in  the  Federal  Government 
for  assuring  maximum  application  of  the  physical  and  engineering  sciences  to 
the  advancement  of  technology  in  industry  and  commerce.  To  this  end  the 
Bureau  conducts  research  and  provides  central  national  services  in  three  broad 
program  areas  and  provides  central  national  services  in  a  fourth.  These  are: 
(1)  basic  measurements  and  standards,  (2)  materials  measurements  and 
standards,  (3)  technological  measurements  and  standards,  and  (4)  transfer  of 
technology. 

The  Bureau  comprises  the  Institute  for  Basic  Standards,  the  Institute  for 
Materials  Research,  the  Institute  for  Applied  Technology,  and  the  Center  for 
Radiation  Research. 

THE  INSTITUTE  FOR  BASIC  STANDARDS  provides  the  central  basis  within  the  United 
States  of  a  complete  and  consistent  system  of  physical  measurement, 
coordinates  that  system  with  the  measurement  systems  of  other  nations,  and 
furnishes  essential  services  leading  to  accurate  and  uniform  physical 
measurements  throughout  the  Nation's  scientific  community,  industry,  and 
commerce.  The  Institute  consists  of  an  Office  of  Standard  Reference  Data  and 
a  group  of  divisions  organized  by  the  following  areas  of  science  and 
engineering: 

Applied  Mathematics — Electricity — Metrology — Mechanics — Heat — Atomic 
Physics — Cryogenics2 — Radio  Physics2 — -Radio  Engineering2 — Astrophysics2  — 
Time  and  Frequency.2 


THE  INSTITUTE  FOR  MATERIALS  RESEARCH  conducts  materials  research  leading  to 
methods,  standards  of  measurement,  and  data  needed  by  industry,  commerce, 
educational  institutions,  and  government.  The  Institute  also  provides 
advisory  and  research  services  to  other  government  agencies.  The  Institute 
consists  of  an  Office  of  Standard  Reference  Materials  and  a  group  of  divisions 
organized  by  the  following  areas  of  materials  research: 


Analytical  Chemistry — Polymers — Metallurgy — Inorganic  Materials — Physi 
Chemistry. 
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THE  INSTITUTE  FOR  APPLIED  TECHNOLOGY  provides  for  the  creation  of  appropriate 
opportunities  for  the  use  and  application  of  technology  within  the  Federal 
Government  and  within  the  civilian  sector  of  American  industry.  The  primary 
functions  of  the  Institute  may  be  broadly  classified  as  programs  relating  to 
technological  measurements  and  standards  and  techniques  for  the  transfer  of 
technology.  The  Institute  consists  of  a  Clearinghouse  for  Scientific  and 
Technical  Information,3  a  Center  for  Computer  Sciences  and  Technology,  and  a 
group  of  technical  divisions  and  offices  organized  by  the  following  fields  of 
technology: 


Building  Research — Electronic  Instrumentation — Technical  Analysis — Product 
Evaluation — Invention  and  Innovation — Weights  and  Measures — Engineering 
Standards — Vehicle  Systems  Research. 

THE  CENTER  FOR  RADIATION  RESEARCH  engages  in  research,  measurement,  and 
application  of  radiation  to  the  solution  of  Bureau  mission  problems  and  the 
problems  of  other  agencies  and  institutions.  The  Center  for  Radiation 
Research  consists  of  the  following  divisions: 

Reactor  Radiation — Linac  Radiation — Applied  Radiation — Nuclear  Radiation. 


headquarters  and  Laboratories  at  Gaithersburg,   Maryland,   unless   otherwise 

noted;  mailing  address  Washington,  D.C.  20234. 

2Located  at  Boulder,  Colorado  80302. 

3Located  at  5285  Port  Royal  Road,  Springfield,  Virginia  22151. 


UNITED  STATES  DEPARTMENT  OF  COMMERCE 
C.  R.  Smith,  Secretary 

NATIONAL  BUREAU  OF  STANDARDS   •   A.  V.  Astin,  Director 


NBS 


TECHNICAL  NOTE  446 

ISSUED  JUNE  1968 


PRECISE:  A  Multiple 
Precision  Version  of  Omnitab 


Alfred  E.  Beam 

Computer  Science  Center 

University  of  Maryland 

College  Park,  Maryland   20742 

and 

Joseph  Hilsenrath 

Office  of  Standard  Reference  Data 

Institute  for  Basic  Standards 

National  Bureau  of  Standards 

Washington,  D.C.    20234 


NBS  Technical  Notes  are  designed  to  supplement  the 
Bureau's  regular  publications  program.  They  provide  a 
means  for  making  available  scientific  data  that  are  of 
transient  or  limited  interest.  Technical  Notes  may  be 
listed  or  referred  to  in  the  open  literature. 


For  sale  by  the  Superintendent  of  Documents,  U.S.  Government  Printing  Office,  Washington,  D.C,  20402 

Price  55  cents. 


ABSTRACT 

This  users  manual  describes  PRECISE  -  a 
completely  assembled  interpretive  program  for  the 
IBM  7090/7094  which  enables  the  user  to  carry  out 
arithmetic  operations  and  function  generation  in 
multiple  precision  (accuracy  to  28  significant 
figures).  PRECISE  operates  as  a  sub-monitor  under 
the  IBSYS  or  DC-IBSYS  monitor  systems.  Appendixes 
describe  how  jobs  are  set  up  to  be  run  under  the 
PRECISE  sub-monitor,  and  how  the  system  may  be 
expanded  to  include  new  subroutines.  The  program, 
which  responds  to  instructions  in  the  form  of 
plain  English  sentences  or  contractions  thereof, 
has  provision  for  handling  numbers  out  of  the 
normal  7090/7094  range.  It  handles  numbers  as 
large  as  10  to  the  10  to  the  9  power.  Other 
features  of  the  program  include:  free-field 
input;  a  work-sheet  of  7,500  cells  (3x2500 
computer  words)  which  can  be  dimensioned  by  the 
user  at  run  time  (75  rows  by  100  columns,  300  rows 
by  25  columns,  (etc.);  solution  of  systems  of 
linear  equations  in  as  many  as  85  unknowns; 
flexible  formatting;  tape  handling  facility;  and 
row  and  column  sums.  A  description  of  the  UOM 
Multiple  Precision  Package  (SHARE  Dist.  No.  3081) 
is  included  as  an  appendix. 


Key  Words:  Elementary  functions,  equation 
solver,  double  precision,  linear  equation  solver, 
magnetic  tape  utility  program,  multiple-precision 
computing,  multiple-precision  programming, 
PRECISE,  triple-precision,  user's  manual. 
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1.   INTRODUCTION 


One  of  the  more  troublesome  problems  that  confront  the 
careful  user  of  modern  computers  is  the  loss  of  significance 
resulting  from  round  off  and  other  computing  pitfalls.  In 
many  calculations  rounding  errors  are  serious  sources  of 
annoyance  --  in  some  they  are  downright  fatal.  While  the 
recent  trend  to  build  computers  with  built-in  hardware  for 
double-precision  operations  is  a  decided  help  in  this 
regard,  the  careful  user  of  these  features  must  still  be  on 
guard.  He  must  be  on  guard  for  possible  flaws  in  the 
hardware,  or  in  the  algorithms  and  even,  unhappily,  for 
errors  in  important  constants  used  by  the  compiler  or  the 
conversion  routines. 

The  problem  has  gotten  worse  recently  as  a  consequence 
of  the  fact  that  many  of  the  third  generation  computers  have 
a  shorter  word  length.  As  a  result,  programs  which 
previously  gave  suitable  answers  in  single  precision  now 
must  be  run  in  double  precision. 

We  are  tempted  to  speculate  that  if  the  cloak  of 
anonymity  were  removed  from  commercial  software  systems  and 
each  subroutine  or  program  or  compiler  segment  were  to  carry 
the  by-line  of  its  author  or  authors,  then,  perhaps,  there 
might  be  some  improvement  in  the  situation.  But,  be  that  as 
it  may,  there  is  a  clear  need  for  some  yardsticks  by  which 
the  accuracy  of  computer  results  can  be  judged.  There  is  a 
need  for  a  system  which  can  deliver  correct  answers  to  a 
reasonably  large  number  of  significant  digits  even  when 
handling  exceedingly  large  or  small  numbers. 

The  release  to  SHARE  in  1963  of  a  multiple  precision 
package  (UOM  MPP  SHARE  DIST.  NO.  3081)  by  Alfred  E.  Beam  was 
a  considerable  boon  to  professional  programmers  using  IBM 
7090-94  computers.  In  spite  of  the  existence  of  the  MPP 
package  and  doubtless  other  similar  packages,  the  problem  of 
carrying  out  calculations  in  multiple  precision  involving 
the  elementary  trigonometric  and  transcendental  functions  is 
still  by  no  means  a  trivial  job.  Nor  is  it  easy  even  today 
to  solve  a  large  system  of  linear  equations  (in  say  85 
unknowns)  and  retain  adequate  accuracy. 
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Last  place  "errors"  are  so  much  a  part  of  even  reliable 
mathematical  tables  as  to  cause  L.  J.  Comrie,  a  well-known 
table  maker,  to  write  a  short  piece  entitled  "What  is  an 
Error"  (MTAC*,  V.2,  1943,  pp  284-286)  in  which  he  explains 
that  when  the  seventh,  eight,  and  ninth  places  in  an  entry 
in  a  mathematical  table  are  4,9,9  or  5,0,0,  it  matters 
little  to  the  man  who  wants  only  seven  places  exactly  what 
the  tenth  or  eleventh  place  is.  Thus,  Comrie  continues 
"...on  more  than  one  occasion  I  have  written  to  our  beloved 
editor  saying  'I  have  found. . .errors  of  less  than  one  unit 
in... tables,  but  am  not  sending  them  to  you,  lest  you  should 
be  tempted  to  publish  them." 

Table  makers  are  quite  willing  to  accept  these  last- 
place  or  end  figure  "errors"  because  of  the  tedium  of 
carrying  out  check  calculations  to  three  or  more  figures 
beyond  those  that  they  normally  carry.  PRECISE  carries  out 
calculations  to  many  figures  as  a  matter  of  course.  Thus, 
there  is  really  no  need  to  tolerate  "end-figure"  errors. 

Soon  after  it  became  clear  that  the  philosophy  behind, 
the  organization  and  implementation  of  the  OMNITAB  general- 
purpose  computing  program  on  the  7094  was  sound  enough  to 
attract  a  wide  audience  of  problem  solvers,  whom  even 
FORTRAN  had  not  reached,  we  turned  our  attention  for  a  time 
to  the  design  of  a  comparable  system  for  more  precise 
calculations  than  were  then  possible  in  single  precision. 
This  system  drew  heavily  on  the  multiple  precision  package 
designed  by  one  of  the  authors  to  spare  professional 
programmers  the  tedium  of  writing  painstaking  instructions 
for  the  computer  to  handle  double  and  triple  precision  and 
out-of-range  arithmetic.  This  report  describes  how  the  MPP 
package  has  been  further  employed  to  provide  nonprogrammers 
with  a  computer  tool  for  very  precise  calculations  without 
the  need  to  resort  to  conventional,  and  in  this  instance, 
very  tedious  programming. 


*  Mathematical  Tables  and  Other  Aids  to  Computation   is   now 
called  "Mathematics  of  Computation." 
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The  PRECISE  program  which  is  discussed  here  was 
designed  to  carry  out  arithmetic  operations  and  function 
generation  often  to  as  many  as  28  significant  digits  and  at 
the  very  least  to  21  figures.  Except  when  instructed  to 
increase  the  ranges,  the  program  normally  handles  numbers  x 
in  the  range  10exp(-76)  to  10exp(76)  and  gives  results  to  28 
significant  figures.  The  program  can  also  handle  numbers 
outside  of  the  above  range.  The  greatest  or  smallest  power 
of  10  can  be  as  high  as  plus  or  minus  one  billion.  In  this 
extreme  case  the  results  are  good  only  to  21  significant 
figures. 

PRECISE,  like  its  predecessor,  OMNI TAB*,  is  designed  to 
provide  a  close  parallel  to  the  modus  operandi  in  carrying 
out  calculations  with  a  desk  calculator  and  a  multi-columned 
(and  multi-lined)  worksheet.  While  the  worksheet  in  OMNITAB 
for  the  7094  was  fixed  at  101  rows  by  46  columns,  the  7500 
cells  (3  x  2500  computer  words)  set  aside  for  the  worksheet 
in  PRECISE  can  be  dimensioned  at  the  start  of  each  problem 
at  run  time. 


*J. 

and  R.   J. 
Statistical 
Standards, 
Washington, 


Hilsenrath,  G.  G.  Ziegler,  C.  G.  Messina,  P.  J.  Wals 
Herbold,  OMNITAB  -  A  Computer  Program  F 
and  Numerical  Analysis,   National   Bureau 

Handbook  101,  U.  S.  Government  Printing  Offic 

D.  C.  20402  (March  4,  1966). 
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2.  TOE  PRECISE  PROGRAM 


PRECISE  is  a  completely  assembled  interpretive  program 
for  the  IBM  7090/7094  which  enables  the  user  to  carry  out 
arithmetic  operations  and  function  generation  -  in  multiple 
precision  (accuracy  to  28  significant  figures) .  PRECISE 
operates  as  a  sub-monitor  under  the  IBSYS  or  DC-IBSYS 
monitor  systems.  Appendix  I  describes  how  jobs  are  set  up 
to  be  run  under  the  PRECISE  sub-monitor.  The  program,  which 
responds  to  instructions  in  the  form  of  plain  English 
sentences  or  contractions  thereof,  has  provision  for  han- 
dling numbers  out  of  the  normal  IBM  7090/7094  range.  It 
handles  numbers  as  large  as  10  to  the  10  to  the  9  power. 
Other  features  of  the  program  include: 

a.  free  field  input 

b.  a  work-space  of  about  7500  cells  which  is  normally  set 
to  simulate  a  work-sheet  of  200  lines  by  37  columns  but 
can  be  rediraensioned  by  the  command  NUMBER  OF  ROWS  IS 
xx,  MAXIMUM  COLS  IS  xx,  or  simply  by  NUMBER  1000,7  for  a 
work-sheet  of  1000  lines  and  7  columns 

c.  A  simultaneous  linear  equation  solver  which  can  handle 
as  many  as  85  equations  with  full  precision 

d.  format  free  output  permitting  printing  and  punching  of 
results  either  in  fixed  or  floating  mode  with  9,  14,  17, 
25,  28  figures  decimals  as  the  case  may  be 

e.  formatting  provision  to  provide  a  wide  variety  of  output 
corresponding  to  the  more  common  practice  in  published 
mathematical  tables 

f.  provision  for  column  headings  of  up  to  18  characters 

g.  provision  for  taking  sums  of  rows  and  of  columns  affords 
a  moderately  efficient  means  for  carrying  out  numerical 
integration  -  with  Gaussian  or  Lagrangian  coefficients 

h.   tape  handling  facilities 

The  over-all  logic  and  command  structure  is  similar  to 
that  of  the  OMNITAB*  program.  There  are  however  the 
following  differences 

a.  the  "work-sheet1'  can  be  dimensioned  (the  product  of  rows 
by  columns  cannot  however  exceed  7500) 

b.  format  statements  permit  more  sophisticated  control 
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c.  the  program  is  normally  set  to  handle  numbers  in  the 
range  10~76  to  1076  but  the  range  can  be  extended  by  the 
command  GREATEST  POSER  OF  TEN  IS  xx,  where  xx  may  be  as 
large  as  109 

d-  the  first  use  of  the  command  SUMMARIZE  provides  a  column 
sum  for  every  column  printed;  the  second  use  turns  the 
provision  off,  etc. 


The  program  was  written  in  FAP  (or  MAP)  and  makes  use   of 

an  arithmetic  package  called  MULTIPLE  PRECISION  PACKAGE   UOM 

MPP   for   the   IBM   7090/7 094,   developed   by  Alfred  E.  Beam 

(distributed  to  SHARE  as  Dist  No  3081).   The  PRECISE  program 

responds  to  commands  such  as: 

a.  GENERATE  ARGS  . 0  (. 0  1 ) 1 .  ( 1. ) 50,   IN  COL  8 

b.  READ  COL  1,2,4,7 

c.  sub  1.  from  col  11  store  in  col  5 

D.  DIVIDE  COL  21  BY  5040.   MULT  BY  COL  7  ADD  TO  COL  8 

e.  FLOATING  WITH  28  SIGNIFICANT  FIGURES 

f.  LOGE  OF  COL  3  MULT  BY  -1.   ADD  TO  COL  5 

g.  NEGEXP  COL  1  STORE  IN  COL  2 
h.  PRINT  COL  1,4,8,16 

The  logic  of  the  command  structure  is  such  that  the 
program  also  interprets  the  above  instructions  correctly 
even  if  all  of  the  words  but  the  first  are  omitted.  Thus 
the  above  instructions  can  also  be  written  as  follows: 

GENERA  .  0  (.01)  1.  (1.)  50.   8 

READ  1,2,4,7 

SUB  1. ,     11,5 

DIVIDE  21,5040. ,  7,  8 

FLOATING  28 

LOGE  3,  -1.,  5 

NEGEXP  1,  2 

PRINT  1  4  8  16 

How  this  is  achieved  will  be  discussed  later.  At  this 
point  it  should  be  observed  that  the  presence  of  a  decimal 
point  in  a  number,  identifies  it  as  a  particular  number. 
The  absence  of  a  decimal  point  identifies  it  as  a  column 
number.  The  input  is  free  field,  the  numbers  need  only  be 
separated  by  one  or  more  spaces,  by  a  comma,  or  any 
non-numeric  character  other  than  +,-,E  and  the  decimal 
point.  See  Figure  1  for  examples  of  the  variety  of  mixed  data  formats 
which  PRECISE  accepts. 
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The  following  arithmetic  and  mathematical  operations   are 
provided- 


ADD 

SUBTRACT 

MULTIPLY 

DIVIDE 

NEGATIVE 

SOB 

MULT 

DIV 

SQRT 

RAISE 

RAISEI 

ABSOLUTE 

LOGE 

LOGTEN 

EXP 

NEGEXP 

SIN 

COS 

TAN 

COT 

SIND 

COSD 

TAND 

COTD 

SINH 

COSH 

TANH 

COTH 

ARCSIN 

ARCCOS 

ARCTAN 

ARCCOT 

ASIND 

ACOSD 

ATAND 

ACOTD 

ASINH 

ACOSH 

ATANH 

ACOTH 

The  arithmetic  operations  and  function  generating 
facilities  of  PRECISE  are  augmented  by  a  flexible  input  and 
a  number  of  output  option  which  permit  one  to  print  results 
to  9,  14,  17,  25,  or  28  figures  in  either  fixed  point  or 
floating  point  without  writing  format  statements.  Where 
more  flexibility  is  needed  than  is  afforded  by  the  built-in 
formats,  the  program  recognizes  and  interprets  format 
statements  providing  considerably  more  flexibility  than  is 
available  in  FORTRAN. 

Extensive  tape  handling  facilities  have  been  built  into 
PRECISE  to  provide  a  back-up  store  when  the  problem  begins 
to  tax  the  work-space  of  7500  cells  (3x7500  computer  words). 
The  system  provides  for  mounting  and  loading  of  tapes;  for 
setting  the  read  and  write  density;  for  positioning  tapes 
and  skipping  files  and  records  (in  both  directions).  These 
tape  manipulation  facilities  permit  one  to  transfer  (PUT) 
data  from  the  core  (work  sheet)  in  binary  form  and  to 
retrieve  (GET)  them  again  in  the  same  run  or  in  subsequent 
runs.  The  program  can  read  from  as  many  as  five  tape  units 
and  punch  out  print  on  an  equal  number.  Input  tapes  can  be 
switched  during  a  run  between  any  of  six  logical  tape  unit 
numbers  (2,  3,  4,  5,  9,  or  10). 
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As   the   majority   of 
designed  to  handle  require 
input   data,   the  program 
tion,  after  its  execution, 
the  instructions  at  the  en 
program   is   applied   to 
involve  hundreds  or  even  t 
ic  listing  of  input  is   no 
suppressed   by  the  instruc 
can  be  restored  again  with 


the   problems   which   PR 

a   moderately   small 
automatically  writes  eac 

on  a  scratch  tape  and  t 
d  of  the  calculation, 
a   problem  in  which  the 
housands  of  cards  and  th 
t   deemed   important,   i 
tion  NOLIST.   The  listin 

the  word  LIST. 


ECISE  was 
amount  of 
h  instruc- 
hen  prints 
Where  the 
input  data 
e  automat- 
t  can  be 
g  of  input 


In  the  normal  mode  of  operation,  each  time  the  PRINT 
instruction  is  used,  a  new  page  is  called.  Also  column 
headings  are  supplied  over  each  of  the  columns  printed. 
There  are  cases,  however,  when  the  user  would  prefer  not  to 
start  a  new  page  and  not  to  have  column  headings.  Provision 
is,  therefore,  made  to  suppress  these  by  DONTPAGE  and 
DONTHEAD.  Here  again  these  provisions  can  be  restored  by 
the  use  of  the  instructions  DOPAGE  and  DOHEAD. 

As  the  use  of  tabular  data  is  enhanced  by  arranging  them 
in  blocks  of  10  lines  separated  by  a  blank  line,  this 
provision  is  built  into  the  program.  Again,  if  this  feature 
is  not  wanted,  it  can  be  suppressed  by  the  instruction 
DONTSPACE  and  reinstated  by  DOSPACE. 


SOURCE  PROGRAM  LISTING 


PRECISE 


TITLE  ELEMENTARY  AND  TRANSCENDENTAL  FUNCTIONS  TO  28  FIGURES 
FLOATING  WITH  28  SIGNIFICANT  FIGURES 


GENERATE  .01 ( .01) . 1 ( . 1 ) 1. ( 1. ) 10. { 10. ) 100. ( 100. ) 1000.  IN   1 

SIN  OF  1  STORE  IN  2 


DIVIDE  1.0  BY  2    AND  STORE  IN  4 
HEADING  ON  1  CARD 


PRINT  1,2 
STOP 


END  OF  SOURCE  PROGRAM  LISTING. 


L 


356  LINES  OUTPUT. 
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3.  PRECISE  VOCABULARY.  LOGICAL  STRUCTURE  AND  OPERATIONS 


+  »  A  COLUMN  NUMBER 

*  *  A  FLOATING  POINT  NUMBER 

$  =  A  FIXED  POINT  CONSTANT ( INTEGER  > 

$$=  A  COLUMN  NUMBER  OR  A  FLOATING  POINT  CONSTANT 

INPUT  OPERATIONS 
*4*********  ****************  *************************************** 

PRECISE  ANY  ADDITIONAL  WORDS  IN  CARD  COLUMNS  13-72 

REMARK  ANY  COMMENT  IN  CARD  COLUMNS  7-72 

TITLE  ANY  INFORMATION  IN  CARD  COLUMNS  7-72 

GENERATE  ARGUMENTS  *(*)»(»)...*(»)»  IN  COLUMN  + 


READ 
STOP 


COLUMNS  +♦+♦•••»  AND  + 

ENDS  A  BATCH  OF  PROGRAMS 


************************* ***************************************** 

OUTPUT  OPERATIONS 

FIXED       FORMAT  WITH  $  FIGURES  TO  RIGHT  OF  DECIMAL  POINT 
FLOATING    FORMAT  WITH  $  SIGNIFICANT  FIGURES 


HEAD        COLUMN  +  WITH  CHARACTERS  IN  COLUMNS 
HEADING     ON  $  CARDS  WHICH  MUST  FOLLOW 


55-72  OF  THIS  CARD 


FORMAT      ON  $  CARDS  WHICH 
PRINT COLUMNS  +>  +  ♦>«.» 


MUST  FOLLOW 

AND  4- 


PUNCH 


COLUMNS  +»+♦ 


,  AND  + 


MODIFICATION  OPERATIONS 

GRFATEST  POWFR  OF  TFN  NEEDED  WILL  BE  $  ( LFSS  TPAN  TEN  TO  NINTH) 

NUMBFR      OF  ROWS  =  $  AND  MAXIMUM  COLUMN  NUMBER  =  $ 

SUMMARIZE , 


INTERNAL 
*********** 


OPERATIONS 
******************************************************* 


ROWSUM      $$♦$$, ...,  +  , ♦♦.. •»  +  ,$$  AND  b I  ORE  IN  COLUMN  4- 

RATSFI RAI<;F  SS  TO  $  POWFR  AND  s  fORF  TN  COLUMN  + 

SOLVF       THE  FOLLOWING  SYSTEM  OF  $  LINEAR  EQUATIONS 

CHANGF THE  SIGN  OF  COI  I.1MN  4- 


ERASF 
ADn 


COLUMNS  +»  +  »•••<  AND  4- 

SS  TO  S%     AND  <;TnRF  fM  mi  IIMN  ± 


SUBTRACT    $$  FROM  $$  AND  STORE  IN  COLUMN  + 
Mlil  TTP1  Y $$  RY  $S  AND  STHRF  tN  TOt  UMN  4- 


DIVIDE      $$  BY  $$  AND  STORE  IN  COLUMN  4- 

RATSF $$  TO  THE  $S  POWER  AND  STORE  IN  COLUMN  4- 

ABSOLUTE    VALUE  OF  $$  AND  STORE  IN  COLUMN  4- 

ACOSD OF  $$  AND  STORE  IN  COLUMN  4-  (ARC  COMNE-*-DEGREES ) 


ACOSH 
ACOTH 


OF  $$  AND  STORE  IN  COLUMN  4- 
OF  $$  AND  STORE  IN  COLUMN  4- 


(HYPERBOLIC  COSINE) 
(ARC  COTANGENT — DEGREES) 


ACOTH 
ARCCOS 


OF  $$ 
OF  $$ 


AND 
AND 


STORE 
STORE 


IN  COLUMN  + 
IN  COLUMN  4- 


(HYPERBOLIC  ARC  COTANGENT) 
(ARC  COSINE — RADIANS) 


ARCCOT      OF  $$  AND  STORE  IN  COLUMN  + 
ARCSTN      OF  $$  AND  STORE  IN  COLUMN  4- 


(ARC  COTANGENT — RADIANS) 
(ARC  SINE — RADIANS) 


ARCTAN      OF  $$  AND  STORE  IN  COLUMN  4-  (ARC 
■ASJND OF  SS  AND  STORF  TN  COLUMN  4-  (ARC 


TANGENT — RADIANS) 
S I  NE^QEGBEE£JL__ 
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ASINH 
ATAND 


OF  $$  AND  STORE  IN  COLUMN  + 
OF  S3  AND  *jr>BF     TN  mi  UMN  + 


HYPERBOLIC  ARC  SINE) 
ARC  TANGENT— DEGREES) 


ATANH 
XOS 


OF    $$    AND    STORE    IN    COLUMN    + 
OF     S«E     AND     <;TORF     TN     COt  1 IMN     + 


HYPERBOLIC    ARC    TANGENT) 
COMNF — RADIANS) 


COSD 
COSH 


OF  $$  AND  STORE  IN  COLUMN  -*• 
OF  $$  AND  STORF  TN  COI  UMN  + 


COSINE — DEGREES) 
HYPFRBGLIC  COSINE) 


COT 
jCQXD- 


OF  $$  AND  STORE  IN  COLUMN  + 
OF  S3  AND  STORF  TN  TOI  UMN  + 


COTANGENT — RADIANS) 
cotangent— DEGREES ) 


COTH 
FXP 


OF  $$  AND  STORE  IN  COLUMN 

fiF     **     AND     STORF     TN     mi  UMN 


HYPERBOLIC  COIANGENI") 
FXPONFNITAI   FUNCIION) 


LOGF        OF  $$  AND  STORE  IN  COLUMN  + 
1  nr,TFN OF  $S  AND  STORF  TN  COLUMN  4- 


NATURAL  LOGARITHM) 
COMMON  LOGARITHM) 


NEGATIVE    OF  $$  AND  STORE  IN  COLUMN  + 
NFftFvP OF  $$  AND  STORF  IN  CO!  UMN  + 


CHANGES  SIGN  AND  STORES) 
NFGATIVE  EXPONENTIAL) 


SIN 
ST  NO 


OF  $$  AND  STORE  IN  COLUMN 
of  ss  ANT)  STOPF  TN  coiiimn 


SINE — RADIANS) 
stnf nFGRFFS) 


SINH 
SORT 


OF  $$  AND  STORE  IN  COLUMN  + 
OF  S*  AND  STORE  IN  COLUMN  4- 


HYPERBOLIC  SINE) 
SQUARE  ROOT) 


TAN 
TAND 


OF  $$  AND  STORE  IN  COLUMN  + 
OF  *$  AND  STORE  IN  COLUMN  4- 


TANGENT — RADIANS) 
TANGENT— DEGRFES) 


TANH 
ADD 


OF  $$  AND  STORE  IN  COLUMN 
$H  TO  **»  MULTIPLY  BY  ti1i» 


4-  (HYPERBOLIC  TANGENT) 
AND  ADD  TO  CQI UMN  + 


SUBTRACT 
MULTIPLY 


$$  FROM  $$,  MULTIPLY  BY  $ 
%%    BY  **»  MULTIPLY  BY  $* » 


$♦  AND  ADD  TO  COLUMN  4- 
AND  ADD  TO  COLUMN  + 


DIVIDE 
RAISF 


$$  BY  $$♦  MULTIPLY  BY  $$♦ 
$$  TO  $$  POWERt  MULTIPIY 


AND  ADD  TC  COLUMN  + 
BY  $1  AND  ADD  TO  COI  I.1MN  + 


ABSOLUTE    VALUE  OF  $$,  MULTIPLY  BY  $$  AND  ADD  TO  COLUMN  + 
ACOSD OF  *%*    Ml.ll  TTPI  Y  BY  SS  AND  ADD  TO  COl  UMN  + 


ACOSH 
ACOTD 


OF  $$ 
OF  S3 


MULT 
MULT 


IPLY 
IPI  Y 


BY  $$ 
BY  %% 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
COt  UMN  •>• 


ACOTH 
ARCCOS 


OF  $$ 
OF  S3 


MULT 
MULT 


IPLY 
IPI  Y 


BY  $$ 
BY  SEE 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
cm  UMN  ± 


ARCCOT 
ARCSTN 


OF  $$ 
OF  ** 


MULT 
MULT 


IPLY 
IPI  Y 


BY  $$ 
RY  ffi* 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
cm  iimn  ± 


ARCTAN 
ASINH 


OF  $$ 
OF  tt<£ 


MULT 
Mill  T 


IPLY 
IPI  Y 


BY  $S 
BY  $S 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
rot  UMN  ± 


ASINH 
ATANn 


OF  $$ 
OP  $$ 


MULT 
Mill  T 


IPLY 
TPLY 


BY  $$ 
BY  $S 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
COLUMN  4- 


ATANH 
.COS 


OF  $$ 
OF  ** 


MULT 
Mm  T 


IPLY 
TPI  Y 


BY  $$ 

RY  tt 


AND 
AND 


ADD  TO 
Ann  to 


COLUMN  + 
COLUMN  4- 


COSD 
COSH 


OF  $$ 
OF  %% 


MULT 
Ml.ll  T 


IPLY 
IPI  Y 


BY  $$ 
BY  $$ 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
COHJMN  + 


COT 

COTD 


OF  $$ 
OF  $% 


MULT 
MUI  T 


IPLY 
IPI  Y 


BY  $$ 
BY  S$ 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
COLUMN  + 


COTH 
FXP 


OF  $$ 


MULT 
Mi  H  T 


IPLY 
TPI  v 


BY  $$ 
BY  $S 


AND 
ANn 


ADD  TO 
Ann  TO 


COLUMN  4- 
COI  UMN  4- 


LOGF 
I.  OATEN 


OF  $$ 
OF  $S 


MULT 
Mill  T 


IPLY 
TPI  Y 


BY  $$ 
RY  SS 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
COLUMN  + 


NFGATIVE 
NFCFvP 


OF  $$ 
OF  $S 


MULT 
Milt  T 


IPLY 
TPLY 


BY  $$ 
BY  $$ 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  4- 
COLUMN  + 


SIN 
SIND 


OF  $$ 
OF  $$ 


MULT 
MULT 


IPLY 
TPLY 


BY  $$ 
BY  $$ 


AND 
AND 


ADD  TO 

ADD  TO 


COLUMN  + 
COLUMN  4- 


SINH 
SORT 


OF  $$ 
OF  $$ 


MULT 
MULT 


IPLY 
IPLY 


BY  $$ 
BY  $$ 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
COLUMN  4- 


TAN 
TAND 


OF  $$ 
OF  $$ 


MULT 
MULT 


IPLY 
TPLY 


BY  $$ 
BY  $$ 


AND 
AND 


ADD  TO 
ADD  TO 


COLUMN  + 
COLUMN  ± 


TANH 


OF  $$ 


MULTIPLY  BY  $$  AND  ADD  TO  COLUMN  + 
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********************************************************************** 

TAPE  MANIPULATION  OPERATIONS($  MAY  BE  LOGICAL  UNIT  2,3»4»9*  OR  1U) 
*************************  *******************  **  ************************ 

REWIND  TAPES  $♦$»•••»$ 

UNLOAD  TAPES  $»$♦.««>$ 

SETLOW  DENSITY  ON  TAPES  $♦$♦♦••»$ 

SETHIGH  DENSITY  ON  TAPES  $♦$»»»»»$ 

FNDFTL^  TAPES  $»$*••••$  ~ 

POSITION  TAPE  S»S  FILES  AND  $  RECORDS  BACKWARDS*  THEN  $  FILES 

AND  $  RECORDS  FORWARDS 
***************************  *****  *******  **  ********  ********************* 

BINARY  TAPE  INPUT/OUTPUT  OPERATIONS(S  MAY  BE  UNIT  2>3»4»9»  OR  10) 
********************************************************************** 

PUT         ON  TAPE  $»  COLUMNS  +,+»...♦+ 


GET 


FROM  TAPE  $♦  COLUMNS  +»+♦...»+ 


********************************■-»-******■*  *****  ***  ********************** 


OTHER  BCD  TAPE  OPERATIONS 


********************************************************************** 
PRNT?       COLUMNS  +♦...♦+  (SAME  AS  PRINT  EXCEPT  TAPE   2  IS  WRITTEN) 

PRNT?       COLUMNS  +»•, 
PRNT4       COLUMNS  +... 

,.♦  +  (SAME  AS  PRINT  EXCEPT  TAPE   3  IS  WRITTEN) 
,••+  (SAME  AS  PRINT  EXCEPT  TAPE   4  IS  WRITTEN) 

PRNTO       COLUMNS  +  ».. 
PRNT10      COLUMNS  +*.. 

,.,+  (SAME  AS  PRINT  EXCEPT  TAPE   9  IS  WRITTEN) 
,..+  (SAME  AS  PRINT  EXCEPT  TAPE  10  IS  WRITTEN) 

PNCH2       COLUMNS  +♦., 
PNCH?       COLUMNS  +..< 

,.♦+  (SAME  AS  PUNCH  EXCEPT  TAPE   2  IS  WRITTEN) 
,.,+  (SAME  AS  PUNCH  EXCEPT  TAPE   3  lb  WRITIEN) 

PNCH4       COLUMNS  +♦.. 
PNCH9       COLUMNS  +».. 

►  .♦  +  (SAME  AS  PUNCH  EXCEPT  TAPE   4  IS  WRITTEN) 
...+  (SAME  AS  PUNCH  EXCEFT  TAPE   9  IS  WRITTEN) 

PNCHjO      COLUMNS  +  »., 
READ?       COLUMNS  +♦.. 

,♦,+  (SAME  AS  PUNCH  EXCEPT  TAPE  10  IS  WRITTEN) 
.  .♦+  (SAME  AS  READ  EXCEPT  DATA  IS  FROM  TAPE   2) 

READ?       COLUMNS  +♦., 
RFADA       COLUMNS  +  .... 

, .♦+  (SAME  AS  READ  EXCEPT  DATA  IS  FROM  TAPE   3) 
..»+  (SAME  AS  RFAD  EXCEPT  DATA  IS  FROM  TAPE   4) 

READO       COLUMNS  +»•< 
READiO      COLUMNS  +  ».. 

..♦+  (SAME  AS  READ  EXCEPT  DATA  IS  FROM  TAPE   9) 
..♦+  (SAMF  AS  RFAD  EXCEPT  DATA  IS  FROM  TAPE  10) 

SWITCH 
TRANSFER 


INPUT  T 
FOLLOWI 


APE  TO  BE 
NG  DATA  TO 


LOGICAL  TAPE  $(MAY  3E  2»3»4»5*9»  OR  10) 
LOGICAL  TAPE  $(MAY  BE  7»3»4»9>  OR  10) 


ENDTRANSFE 
*********** 


R  OF  DAT 
******** 


AtTERMINAT 
********** 


ES  THE  TRANSFER  OPERATION) 
***************************************** 


OUTPUT  CONTROL  OPERATIONS 


******************* ************#»***********^W**WW***W ******* ********* 


NOLI  ST 


SUPPRESSES  LISTING  OF  THE  SOURCE  PROGRAM 


LIST  RESTORES  LISTING  OF  THE  SOURCE  PROGRAM 

DONTHEAD  SUPPRESSES  PAGE  HEADING  DURING  THE  PRINT  OPERATION 

DOHFAD  RESTORES  PAGE  HEADING  DURING  THE  PRINT  OPERATION 

DONTPAGE  SUPPRESSES  BEGINNING  A  NEW  PAGE 


DOPAGE 
DONT^PACE 


RESTORES  BEGINNING  OF  NEW  PAGES 

SUPPRESSES  BLANK  LINES  BETWEEN  TEN  LINE  BLOCKS 


DOSPACE     RESTORES  BLANK  LINES  BETWEEN  TEN  LINE  BLOCKS 
********************************************************************** 
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4.     DESCRIPTION  OF  OPERATIONS   (GENERAL) 

PRECISE:  This      card      must    be   first    in   each    program.      Card 

columns  13-72  are  placed  in  the  page  line  which 
heads  each  page  unless  it  is  overwritten  by  use 
of    the    TITLE    card. 


REMARK: 


Card  columns  7-72  are  printed  on  the  output  tape 
at  the  time  this  card  is  read. 


TITLE: 


GENERATE: 


READ: 


Card  columns  7-72  of  this  card  are  inserted  into 
the  page  line. 

A  column  is  generated  as  specified  and  the 
number  of  rows  is  reset  to  be  the  number  of 
arguments  generated. 

Data  is  read  into  specified  columns,  one  row  per 
card  and  the  number  of  rows  is  reset  to  be  the 
number  of  data  cards  read. 


FIXED: 


FLOATING: 


STOP:  This  card  must  be  the  last  of  a  set  of  PRECISE 
programs  and  is  the  signal  to  send  control  to 
the  monitor  system. 

The  print  format  is  set  to  be  fixed  and  the 
number  of  decimals  to  the  right  of  the  decimal 
point  is  set  to  be  9,  14,  17,  25,  or  28,  as 
specified. 

The  print  format  is  set  to  be  floating  and  the 
number  of  decimals  to  the  right  of  the  decimal 
point  is  set  to  be  9,  14,  17,  25,  or  28  as 
specified. 

Table  A  summarizes  the  output  options  which  are  available  by 
the  use  of  the  FIXED  and  FLOATING  commands.  The  first 
format  in  the  table  is  the  one  which  is  set  initially,  and 
it  remains  in  force  until  a  FIXED,  FLOATING,  or  FORMAT 
command  is  encountered.  APPENDIX  II  gives  a  general  des- 
cription of  formats. 


HEAD: 


An  18  character  heading  for  the  specified  Column 
is  taken  from  card  columns  55-72  of  the  HEAD 
card.  Provision  was  made  to  head  only  columns 
1-63  so  one  should  not  print  a  column  whose 
number  is  greater  than  63  unless  the  HEADING 
card  has  been  used. 
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HEADING 


The  heading  line  is  blanked  and  the  first  72 
characters  are  replaced  by  card  columns  1-72  of 
the  first  card  following  the  HEADING.  If  2 
heading  cards  are  specified  then  heading  line 
characters  73-120  are  replaced  by  card  columns 
1-48  of  the  second  card  following  the  HEADING 
card. 


FORMAT: 


The  format  area  is  filled  with  card  columns  1-72 
of  the  specified  (1  to  4)  number  of  cards  which 
follow  the  FORMAT  card.  Blanks  are  not  allowed 
except  in  Hollerith  fields  and  the  format  must 
begin  in  card  column  1  of  the  first  card.  The 
first  blank  not  in  a  Hollerith  field  terminates 
the  format.  A  format  description  is  given  in 
APPENDIX  II. 


PRINT: 


PUNCH 


GREATEST; 


NUMBER; 


The  specified  columns  are  written  on  the  regular 
output  tape  according  to  the  current  format. 

The  specified  columns  are  written  on  the  regular 
punch  tape  according  to  the  current  format. 

The  parameter  $  =  P  causes  the  entire  package 
to  be  set  to  handle  numbers  (N)  in  the  range 
10**-P<) N | <10**P  If  this  operation  is  not  speci- 
fied then  P  =  76.  If  given,  it  must  be  given 
before  any  computation,  column  input  or  output. 
P  can  be  increased  by  a  following  program  but  P 
cannot  be  decreased  within  the  same  job.  In  the 
initial  condition  when  p=76,  computation  is 
carried  out  using  28  significant  decimal  digits. 
When  p  is  increased,  the  accuracy  is  decreased. 
The  approximate  number  of  significant  decimal 
digits  in  effect  is  30-Q  where  Q  is  the  closest 
power  of  ten  for  p.  Since  p  may  be  set  to  any 
value  in  the  range  from  76  to  1  billion,  the 
accuracy  of  calculations  will  be  from  28  to  21 
significant  decimal  digits. 

The  row  and  column  parameters  can  be  set  at  any 
time  but  the  product  of  the  two  dimensions  must 
not  exceed  the  matrix  storage  capacity. 


SUMMARIZE 


The  column  summing  -switch  is  turned  on  if  off, 
or  off  if  on.  If  the  summing  switch  is  on  when 
a  PRINT  is  executed  then  the  sum  of  the  printed 
columns  will  also  be  printed. 
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5.  DESCRIPTION  OF  THE  INTERNAL  OPERATIONS 

For  description  of  the  internal  operations,  use   is   made 
of  the  following  definitions: 

1)  k  is  the  current  column  length. 

2)  n  is  the  number  of  arguments  specified  on  a  card. 

3)  C (i)  is  the  ith  argument,  i  =  1,2,. ..,n. 

**)   Cc(i)  ]  is  tne  column  of  numbers  specified  by  C  (i) 

If  C(i)  is  not  a  column  number,  [ C (i)  ]  can  be  thought 
of  as  a  column  of  numbers  with  all  elements  equal  to 
C(i) 

5)  Algebraic   symbols   are   used   to  indicate  operations 
repeated  for  each  element  of  the  column. 

6)  [ H  ]< — £2]  means  "column  »  is  replaced  by  column  2". 

The  table  below  gives  an  operation  name   followed   by   a 
description  of  the  function  of  that  particular  operation. 


GHOOP  I 

fiOHSUH    [C(n)  ]< — £C(1)  ]♦ +[C(n~1)  j 

EAISEI   £C(n)  ]< — £C(1)  ]**C{2)  ,  C(2)  must  be  an  integer 

CHANGE    [C(2)]< — [C{1)] 

ERASE    £C(j)  ]< — 0,  j  =  1,2, ...,n 

SOLVE    This   operation   must   be   followed  by  C(1){C(1)+1) 

numbers  punched  (starting  in  column  2)  on   as   many 

cards   as   necessary.    Each  card   must   have   the 

character  D  punched  in  column  1.    Numbers   are  of 

floating   point   type   and  are  separated  by  commas. 

The  first  blank  encountered   terminates   the   card. 

An   example   of  the  operation  and  data  for  solution 

of  the  two  linear  equations 

-3x  ♦  9y  =  -6 

-5x  -  y  =  4 

is  as  follows: 

SOLVE  2  EQUATIONS 

D-3.E0,9.0,-6 

D-00.05E2,-1.0 

D.00004E+5 

All  of  the  above  data  could  have   been   punched   on 

one  card.   The  matrix  and  solutions  are  printed  and 

then   all   columns  are  cleared  before  proceeding  to 

the  next  operation. 


13- 


GBOOP  II  (BINARY  OPERATORS)  * 


ADD 

SUBTRACT 

MULTIPLY 

DIVIDE 

RAISE 


[C(n)  ]<—  £C(1)  J+[C(2)  ] 
£C(n)  ]<— [C(2)  ]-[C(1)  ] 
£C(n)  ]<— [C(1)  ]*£C(2)  ] 
£C(n)  ]<— [C(1)  ]/[C(2)  ] 
[C(n)  ]<--[C(1)  ]**[C{2)  ] 


The  5  operations  above  require  that  the  number  of 
arguments  n  be  3.  There  is  a  corresponding  operation  for 
each  of  these  when  an  additional  argument  is  specified. 

Let  [FJ  be  the  resultant  column  for  any  of  the  above 
operations.  Then  the  column  result  [C(n)]  for  an  additional 
argument  can  be  described  as 

[C(n)  ]<-- £C(n)  ]+£F]*£C(n-1)  ] 

♦Note:  DIV  may  be  used  rather  than  DIVIDE,  MULT  or  MPY  may 
be  used  rather  than  MULTIPLY,  and  SUB  may  be  used 
rather  than  SUBTRACT. 
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GROUP  III  (UNARY  OPERATORS) 


ABSOLUTE 

NEGATIVE 

EXP 

NEGEXP 

LOGE 

LOGTEN 

SQRT 

SIN 

SIND 

COS 

COSD 

TAN 

TAND 

COT 

COTD 

ARCTAN 

ATAND 

ARCCOS 

ACOSD 

ARCSIN 

ASIND 

ARCCOT 

ACOTD 

SINH 

COSH 

TANH 

COTH 

ASINH 

ACOSH 

ATANH 

ACOTH 


[c 

[n)  ] 

C  ci 

n)  ] 

[  Cj 

n)  1 

[C 

[n)  ] 

[  CI 

[n)  ] 

[C 

n)  ] 

[Cj 

n)  ] 

I  Cl 

[n)  ] 

[C  I 

n)  ] 

[C 

;n)  ] 

E c 

[n)  3 

[Cl 

n)  3 

[  Ci 

n)  ] 

[CI 

n)  ] 

[  c 

n)  3 

L  c\ 

n)  3 

I  c 

n)  3 

[C 

!n)  ] 

[Cl 

;n)  3 

[  Cj 

;n)  3 

[Ci 

;n)  3 

[Ci 

;n)  3 

[C 

n)  3 

[C| 

[n)  3 

[C 

|n)  3 

[C 

n)  ] 

[  c 

n)  3 

[  c 

n)  j 

[C 

n)  3 

[c 

(n)  3 

[C 

[n)  3 

<- 

-![C(1)  31 

<- 

-  -[C(1)  3 

<- 

-e**[C{1)  J 

<- 

-e**-[C(1)  3 

<- 

-log  [C(1)  3 

<- 

-common  logarithm 

<- 

-the  square  root  o 

<- 

-sin[C(1)  3/  radian 

<- 

-sin[C(1)  ]r  degree 

<- 

-cos[C(1)  3/  radian 

<- 

-cos[C(1)  3#  degree 

<- 

-sin[C(1)  3/cos[C{1 

<- 

-sin[C(1)  3/cos[C{1 

<- 

-cos[C(1)  3/sin[C  (1 

<- 

-cos[C(1)  3/sin[C  (1 

<- 

-tan-i[C{1)  3, 

resu 

<- 

-tan-*[C(1)  3, 

resu 

<- 

-cos~1[  C  ( 1)  3f 

resu 

<- 

-cos~l[  C  ( 1)  ], 

resu 

<- 

-sin-i[C{1)  3, 

resu 

<- 

-sin-l[C{1)  3, 

resu 

<- 

-cot~i[C  (1)  3, 

resu 

<- 

-cot~*[  C(1)  3, 

resu 

<- 

-sinh[C  (1)  3 

<- 

-cosh[C  (1)3 

<- 

-tanh[C  (1)  j 

<- 

-coth[C  (1)  3  = 

1/ta 

<- 

-sinh-i[C  (1)  3 

<- 

-cosh-i[C  {1) 3 

<- 

-tanh-»[C  (1)  3 

<- 

-coth~i[C  (1)  3 

of  [C(1)  3 
f  [C(i)  3 

argument 

argument 

argument 

argument 
)  3r  radian  argument 
)  3#  degree  argument 
)  3#  radian  argument 
)  J,  degree  argument 
It  in  radians 
It  in  degrees 
It  in  radians 
It  in  degrees 
It  in  radians 
It  in  degrees 
It  in  radians 
It  in  degrees 


nh[C(1)  3 


The  above  operations  (Group  III)  require  that  the  number 
of  arguments  n  be  2.  There  is  a  corresponding  operation  for 
each  of  these  when  an  additional  argument  is  specified. 

Let  [F3  be  the  resultant  column  for  any  of  the  above 
operations.  Then  the  column  result  [C(n)j  for  an  additional 
argument  can  be  described  as 

[C(n)  ]<— [C(n)  3+[F]*[C(n-1)  3 
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6.  BINARY  TAPE  OPERATIONS  AND  TAPE  MANIPULATION 


The  following  operations  are  designed  for  positioning  of 
tapes,  writing  tape  marks,  saving  of  temporary  results,  and 
for  reading  of  data  which  was  previously  written  on  tape. 
In  the  descriptions  T  means  a  logical  tape  number  and  T  may 
have  integral  values  equal  to  2,3,4,9,  or  10. 

REWIND  TAPES  T(1),T(2),... 

This  operation  causes  all  tapes  specified  to  be  rewound. 

UNLOAD  TAPES  T  (1)  ,T  (2)_£.  ,. 

This  operation  causes  all  tapes  specified  to  be  rewound  and 
unloaded. 

SETLOtf  DENSITY  ON  TAPES  T  f  1)  ,T  <2)  .  .  .  . 

This  operation  causes  all  tapes  specified  to  be  set  to  low 
density 

SETHIGH  DENSITY  ON  TAPES  Tf1),T(2). 

This  operation  causes  all  tapes  specified  to  be  set  to  high 
density. 

ENDFILE  TAPES  T  (  1)  ,T  (2)  ,  ■  .  . 

This  operation  causes  a  tape  mark  {end  of  file)  to  be 
written  on  all  specified  tapes. 

POSITION  TAPE  Tfa.b,c.d 

This  operation  causes  tape  T  to  be  positioned  as  follows: 

1)  Tape  T  is  moved  backwards  over  a  tape  marks. 

2)  Tape  T  is  moved  backwards  over  b  records. 

3)  Tape  T  is  moved  forward  over  c  tape  marks. 

4)  Tape  T  is  moved  forward  over  d  records. 

Any  of  the  numbers  a,b,c,  or  d  may  be  zero,  but  may  not  be 
omitted  from  the  operation. 

PUT_0N  TAPE  T,.  COLUMNS  C(1).C(2), ,C(n) 

This  operation  causes  all  specified  columns  Ci  (1  record  per 
column)  to  be  written  on  tape  T.  If  no  Ci  are  specified 
then  the  entire  working  area  is  written  as  one  record  on 
tape  T. 
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GET  FROM  TAPE  T,  COLUMNS  C(1),C(2), ,C  (n) 

This  operation  causes  all  specified  columns  (1   record   per 

column)   to  be  replaced  by  columns  obtained  from  tape  T-   If 

no  C±   are  specified  then  the  entire  working  area  is  replaced 
by  the  next  record  from  tape  T.   All  information   read   via 

the  VGET*  operation  should  have  been  put  on  the  tape  via  the 
•POT'  OPERATION. 


OTHER  BCD  TAPE  OPERATIONS 

The  following  operations  are  mainly  useful  for  special 
purpose  jobs,  such  as  building  or  reading  special  input 
tapes  or  preparing  special  output  tapes.  In  the  description 
of  the  operation  T  is  an  integral  logical  tape  number  and 
the  only  legal  values  for  T  are  2,3,4,9,  or  10.  Operation 
names  terminating  with  i  means  that  the  operation  name  is 
defined  for  i  =  2,3,4,9,  or  10  depending  on  the  value  of   i. 

PRNTi     COLUMNS  C  (1)  fC(2)  ,..,,C(n) 

This  operation  is  the  same  as  'PRINT*  except  the  specified 
columns  are  printed  on  logical  tape  i  rather  than  the 
regular  output  tape, 

PNCHi     COLUMNS  C  ( 1)  ,  C  (2)  f ,  .  .  f  C  (n) 

This  operation  is  the  same  as  'PUNCH*  except  the  specified 
columns  are  printed  on  logical  tape  i  rather  than  the 
regular  punch  tape. 

READi     COLUMNS  C  (1j  gC{2)_g.  .  ,  #-CJn)_ 

This  operation  is  the  same  as  'READ'  except  the  specified 
columns  are  read  from  tape  i  rather  than  the  regular  input 
tape. 

SWITCH  INPUT  TAPE  TO  BE  T 

This  operation  causes  the  next  operation  or  data  card  (and 
following  cards)  to  be  read  from  tape  T  rather  than  the 
regular  input  tape. 

TRANSFER  NEXT  STATEMENTS  TO_TAPE .T 

This  operation  causes  all  following  cards  to  be  transfered 
to  tape  T  without  any  attempt  to  interpret  the  cards.  All 
cards  up  to,  but  not  including,  a  *ENDTRA*  operation  will  be 
transfered. 
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ENDTRANSFER  OF  INSTRUCTIONS 

This   operation   terminates  the  transfer  of  cards  to  another 
tape  which  was  initiated  by  the  'TRANSFER'  operation. 

OUTPUT  CONTROL  OPERATIONS 

The  following  operations  give  the  user  some  control   over 
the  heading,  spacing,  and  program  listing. 

NOLIST 

This  operation   suppresses   listing   of   the  PRECISE  source 
program. 

LIST 

This  operation  restores  the  listing  of   the   PRECISE   source 
program.   LIST  is  initially  the  mode  of  listing. 

DONTHEAD 

This  operation   suppresses  heading  of  a  page  when  the  PRINT 
operation  is  encountered. 

DOHEAD 

This  operation  restores   heading   of   pages-    This   is   the 
initial  mode. 

DONTPAGE 

This  operation   suppresses  the  beginning  of  a  new  page  when 
the  PRINT  operation  is  encountered. 

DOPAGE 

This  operation  restores  the  paging  of   output,   and   is   the 
initial  mode. 

DONTSPACE 

This   operation   suppresses  the   printing   of   a  blank  line 
between  each  block  of  ten  lines. 

DQSPACE 

This  operation  restores  the  printing  of  a  blank  line  between 
each  block  of  ten  lines.   This  is  the  initial  mode. 
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7.  RESTRICTIONS  ON  THE  OPERATIONS 


1)  The  operation  word  must  be  the  first  character  punched 
on  a  card. 

2)  Except  for  the  PRECISE,  REMARK,  HEAD,  TITLE,  FORMAT 
data,  and  D-type  data  cards;  spacing  is  arbitrary, 
anywhere  within  the  72  columns  of  a  card.  PRECISE, 
REMARK,  HEAD,  and  TITLE  should  be  punched  starting  in 
card  column  1. 

3)  All  operation  words  must  be  at  least  6  characters  long 
or  separated  from  following  information  by  a  blank. 

4)  An  argument  Ci  is  composed  of  a  continuous  string  from 
the  set  of  characters  0,1,. ..,9,  decimal  point,  E,  ♦  ,  or 
-.  Any  character  (s)  o';her  than  this  set  acts  as  an 
argument  separator.  Ci  will  be  a  multiple  precision 
number  if  either  of  the  characters  E  or  ._  appear  in  the 
field.  Ci  will  be  an  integer  if  neither  E  nor  ._ 
appear. 

5)  In  most  operations  an  integer  argument  must  be  less  than 
or  equal  to  the  maximum  column  number  allowable.  Excep- 
tions to  this  rule  are  GENERATE,  GREATEST,  NUMBER, 
FORMAT,  FIXED,  FLOATING,  HEAD,  HEADING,  SUMMARIZE, 
SOLVE,  REMARK,  TITLE,  RAISEI,  REWIND,  UNLOAD,  SETLOW, 
SETHIGH,  ENDFILE,  PUT,  GET,  and  POSITION. 
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8.  ERRORS  AND  ERROR  DETECTION 


The  following  are  considered  errors  by  PRECISE,  and  cause 
execution  of  a  program  to  be  terminated  whenever  such  an 
error  is  detected. 

1)  An  illegal  operation  word. 

2)  A  colum  number  which  is  too  large. 

3)  Illegal   data. 

4)  Overflow  during  computation  or  input  conversion.  Over- 
flow occurs  when  a  number  is  computed  which  is  too  large 
to  be  represented  within  the  currently  defined  binary 
number  format. 

5)  Division  by  zero. 

6)  Output  FORMAT  errors. 

Whenever  one  of  the  above. errors  is  detected,  one  or  more 
messages  are  printed,  and  execution  of  the  rest  of  the 
PRECISE  commands  is  terminated;  however,  the  remainder  of 
the  commands  are  scanned  to  determine  legality  of  the 
operation  words. 


21 


1 — 

CARD  XXXXXXX 
THIS  POINT 

1— 1 

o 
o  o 

O   LU 

LU  *-* 
•4-  ^ 

1—1 

o 
o 
o 
o 
o 
o 
o 
o 
o 
o 
o 

o 

i— 4 

UJ 

o 
o 
o 
o 
o 
o 
o 
o 
o 

• 

0^ 
0> 

a* 

0> 
0* 

Csl 

Csl   >- 
Csl    _J 
Csl    _J 
CJ    < 

Csl  O 

Csl   t-4 
Csl   h- 
Csl    < 

cj  s: 

csi  a 
eg  »— 

X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 

1- 

k— 4 

i— 

eg 

z 

X 

■ 

1 

LU       ' 
2      1 
O       j 

1—      J 

00      ' 
OJ      ] 

h- 

t—i 

CJ  r-« 

o 

o 

o> 

eg  o 

X 

s: 

s: 

>— » 

3 

>~      11 

*£   < 

LU 

o  -* 

o  o 

o 

O 

eg  < 

X 

H— 1 

»— • 

2: 

i- 

-J 

a: 

LO 

o  ~* 

o  o 

o 

o 

eg 

X 

-J 

_i 

»— 1 

(—4 

a 

u-     I 

<  Q 

CO 

o  -4 

o  o 

o 

o 

eg  oo 

X 

_i 

2: 

0 

0    * 

X   LU 

f\J 

•  1— » 

o  o 

o 

0> 

eg 

X 

z 

z 

>— < 

--4 

ro 

ip 

lu  cc 

o 

>J-  ~i 

1  -1 

o  o 

o 

v> 

eg  uj 

X 

s: 

2! 

z 

_j 

0 

Q     1 

LU 

or  Q 

>i-  o 

1     (Ni 

r-4 

•  o 

o 

c* 

eg  z 

X 

x> 

3 

2: 

z 

1— 

Z 

z 

z    1 

H- 

z 

(\J  o 

LU  LU 

• 

o  o 

o 

o> 

eg  >— i 

X 

_j 

_J 

z> 

z 

X 

z 

X 

X 

LU      1 

< 

II    3 

>o  o 

-*  o 

r- 

o  o 

o 

0> 

eg  z 

X 

o 

O 

-J 

X 

3 

>—i 

3 

3 

CC 

X 

•      • 

<\J  LA 

cg 

o  o 

o 

c* 

eg 

X 

o 

o 

0 

3 

_J 

_J 

_l 

LU 

+ 

CO   CO 

•  l—t 

LO 

o  o 

o 

o> 

eg  x 

0 

_J 

O 

O 

a 

*       , 

z 

1   o 

r~ 

r-  >t 

o  o 

o 

0> 

eg  I— 

a 

LU 

O 

O 

O 

0 

*       1 

LU 

*fi>  3: 

i— i 

1  r-t 

o  o 

o 

o 

(NJ     l-H 

a: 

LU 

a: 

LU 

O 

• 

*       1 

o 

^  h- 

• 

o 

r-4 

O" 

eg  :s 

< 

a: 

3 

> 

z 

O 

0 

O 

a 

*       ' 

>— 

ro 

• 

o 

X 

o 

•—4 

0 

X 

1— 

>*- 

1- 

H 

*       | 

Q 

—   00 

t 

• 

* 

o 

t— 

u. 

UL 

1—4 

z 

Z 

Z 

*       1 

z 

*      H- < 

r-4  C\J 

ro  >fr 

ie\  o 

r-  co 

&• 

o 

o  z 

LU 

or 

00 

1—4 

• 

t— 4 

1—4 

*       1 

< 

INI 

>  00 

I 

1 

1 

• 

•  1— • 

1— 

00 

i— • 

LU 
O 

1— 1 

o 

X    jE 

o 

< 

a 

O 

• 

v» 

3:  a 

o  o 

o  o 

o  o 

o  o 

o 

o 

sO  o 

UJ 

1-^ 

0 

M 

1 

1 

>  or 

3 

o  o 
o  o 
o  o 

o  o 
o  o 
o  o 

o  o 
o  o 
o  o 

o  o 
o  o 
o  o 

o 
© 
o 

o 
o 
o 

m  ^j 

1     UL 

UJ 

a. 

z 

1—4 

O 

?E 

ro 

• 

l-t 

00  z 

o  o 

o  o 

o  o 

o  o 

o 

o 

r-*    00 

UJ 

< 

O 

~-> 

*c 

or 

o  o 

o  o 

o  o 

o  o 

o 

o 

eg  •-• 

X 

O 

z 

z 

r-4 

o 

a  x 

o  o 

o  o 

o  o 

o  o 

o 

o 

ro 

H- 

< 

2: 

1— 1 

t 

• 

LU 

a.  3 

o  o 

o  o 

o  o 

o  o 

o 

o 

>t  1- 

3 

3E 

O 

>g- 

X 

o  x 

o  o 

o  o 

o  o 

o  o 

o 

o 

tn  < 

Z 

vC 

j 

2: 

t-4 

o 

z  •-» 

o  o 

o  o 

o  o 

o  o 

o 

o 

o  s. 

o 

0 

3 

O 

— 1 

x  x 

o  o 

o  o 

o  o 

o  o 

o 

o 

r-  a: 

1— 

0 

00 

» 

• 

* 

-j  < 

o  o 

o  o 

o  o 

o  o 

o 

o 

CO  o 

00 

< 

t-t 

O 

*:  s: 

o  o 

o  o 

o  o 

o  o 

o 

0> 

0^  U- 

< 

LU  O 

00 

t 

<*■ 

• 

x 

"3 

o  o 

o  o 

o  o 

o  o 

CO 

t-4 

o 

jC 

a. 

^- 

h- 

O 

* 

O 

r-4 

< 

-H     >* 

o  o 

o  o 

o  o 

o  r- 

1—1 

o  >t 

LU 

'JJ 

>* 

LTl 

• 

LU 

X 

o  o 

o  o 

o  o 

■o  rM 

o 

H-  O 

o 

eg 

o 

O 

cC 

f-« 

00 

»-» 

— «. 

r-4 

r-4 

co 

o  o 

o  o 

o  o 

O  LO 

— 4 

o  o 

<t  -• 

«-« 

^-4 

•-^ 

4» 

0 

• 

• 

1 

• 

1 

u.  z 

o  o 

o  o 

^  r-4 

o  z 

x  u. 

u. 

O 

u. 

u. 

a 

-• 

Q 

•-* 

F-t 

f-4 

•—» 

X 

LU   < 

o  o 

O  ro 

—4 

o  < 

1- 

~-4 

z 

z 

4ta 

CC 

u. 

- 

t 

r-4 

O      j 

r- 

Q 

o  o 

oo  .-i 

o    » 

o 

o 

< 

o 

0 

< 

eg 

< 

0 

O 

• 

O 

* 

<t 

o    - 

o  -> 

r-4 

O  ro 

O  — 

Q 

— 1 

»— » 

1— « 

•• 

O 

t— < 

O 

•4- 

O      1 

C£ 

co  ro 

O    r-4 

o    » 

U.   00 

Z 

00 

«* 

oo 

ro 

00 

ro 

ro 

00 

00 

00 

00 

z    « 

z 

<    •» 

—4 

•  eg 

z 

< 

r-l 

•» 

•* 

•» 

00 

gD  t- 

00 

k- 

h- 

<    1 

LU 

o>  eg 

1     » 

<— <   CO 

LO 

m 

o 

eg 

>* 

C^i 

>J- 

»— 1 

cr 

*  z 

1— 

Z 

Z 

» 

00 

CO     •> 

—t  —4 

o 

i—4 

eg 

^r 

eg 

<t 

eg 

4* 

r-4 

«h 

X 

r-4   LU 

Z 

LU 

UJ 

m    1 

_J 

r»-  »-• 

o  r-4 

Csl 

1 

>r 

00 

^-1 

r-J 

00 

(— 

•«  X 

LU 

X 

X 

*»    | 

*H 

go 

o  o 

— *  ro 

LU 

o 

LCt   00 

uj  x 

00 

X 

00 

-y 

« 

X 

*» 

Z 

X 

>J-  3 

X 

3 

3 

>t  j 

X 

in  oo 

o  o 

o  <-* 

*■ 

+ 

1    z 

U  H- 

z 

1— 

z 

H- 

>i- 

K 

>3- 

2: 

t— 

•  O 

3 

O 

O 

•  1 

-4-  z 

o  o 

o  o 

r-t    LT» 

LL 

o 

LU  2: 

<   ►— I 

2: 

1—4 

X. 

1—4 

►—4 

z> 

►— 4 

ro  3; 

O 

CC 

a: 

ro    J 

ro  x 

o  o 

o  o 

o  -* 

o 

o 

^  3 

-J  3: 

=>  -x 

3 

3: 

00 

3E 

00 

-J 

3: 

•  < 

or 

< 

< 

• 

eg  3 

o  o 

o  o 

o  o 

*-•  r- 

o 

•  -J 

a 

_l 

_l 

_ I 

-J 

a 

UJ 

eg 

< 

eg 

«-t  _j 

o  o 

o  o 

o  o 

O   fH 

CO 

LU 

f*    O 

uj  o 

a 

o 

o 

o 

O 

0 

O 

u 

NJ 

0 

•>  LU 

LU 

LU 

m      1 

O  O 

o  o 

o  o 

o  o 

o  o 

1-4 

c* 

eg  o 

a:  z 

U  2 

u 

z 

o 

Z 

O 

>— 4 

z 

r-4    1- 

LU 

K 

H- 

*~4        | 

UJ 

<-> 

o  o 

o  o 

o  o 

o  o 

o 

r-4 

ro 

1—4 

(—4 

H^ 

>— • 

2: 

a: 

1— 1 

< 

r- 

< 

< 

oo 

*: 

o  o 

o  o 

o  o 

o  o 

o 

o 

>t  K 

UJ  V- 

H- 

1— 

t- 

H- 

K- 

h- 

l~ 

3 

< 

H- 

I-    C£ 

< 

ct 

<X 

t—        | 

>— ( 

a:  o 

o  o 

o  o 

o  o 

o  o 

o 

o 

LO    Z 

-J  < 

Z  < 

z 

< 

z 

< 

Z 

00 

x 

< 

Z   LU 

ac 

UJ 

LU 

Z     1 

o 

<  < 

•     • 

•     • 

•     • 

•     • 

• 

• 

g3  -. 

K-  o 

i— • 

o 

1—1 

a 

►—4 

O 

1— « 

3: 

X 

O 

-<  z 

LU 

z 

z 

•—4          1 

LU 

3L   LU 

1 

1 

I 

1 

1 

P~  ct 

1-4    _J 

QC 

-J 

(X 

-j 

a: 

_J 

cc 

O 

3  -j 

QL   LU 

Z 

UJ 

UJ 

<£    1 

or 

LU   QC 

co  a 

1-  u. 

a. 

LU 

Q. 

u. 

a 

u. 

a. 

or 

to 

u. 

a  0 

LU 

0 

O 

a    J 

a. 

ae 

! 

O 

22- 


Q 

t— 
O 

o 


Z3 

cc 


co 

LU 
CC 

i— i 

u. 
in 

a 


00 

z 
a 


cj 

Z 

Z) 

a. 


2: 

LU 

a 
z 

OJ 

cj 
*/> 

z 
< 


CJ 

*£ 
«1 

>- 
cc 

«. 

z 

LU 

X. 

LU 
-J 

LU   LU 

CO 


a 
o 


• 

o 
o 
o 


o 

a 


o 


o 


*-  cj 
•  o 


a 

LU  LU 
CC  CC 
LU 

H  < 
Z 

13 
O 
CJ 

z 

LU 


a 
0 

X 

3c  a 
o  0 


1— 


O 

z 

< 


X 

z 


X 
co 
a 
o 


Z  3 

LU  C 

LU  I 

CD  H- 

CO  Ou 

<  z 

X  «-• 

z  cj 

a  z 

1-4  M 

H  U 

CJ  3 

Z  -I 

z>  0 

u.  Z 

H 

a  1 

cu  1 

Z  LU 


u. 

LU 

a 


• 
in  »-* 
•r  *~ 
t  ^* 
ir\  • 
cn  — 

LU  •■»• 
f-t      • 

qc-h 
a  wo 
<i  1-1    • 

m  «-« 


CO 


z 

LU 


u 

CC 
CU 


u 

a: 
<i 
cj 


a  — 
a.  t- 


o  o> 

r-  CNJ 

o.  u_ 

cc 
o  u 


u 

a.      z 

eg  cJ 

CO  *-«  CJ 

cu       z 

r—  JL  »-h 

<  Z  Q 

LU  »-4  < 

CC  CO  CU 

o<  I 


CJ 

< 
CJ 


X  NINO  X 

»     •» 

1— 1  ^  CJ 

z 

t-  X  •-• 
Z  CO  Q 

Ct    CJ  CU 

a.  <  x 


LU 
X 

>• 

< 


«L   CU 
CJ 

CC  <l 
CJ   -I 

a. 


_j  ic 

LL.  <L 

CC  r- 
LU 

->  > 

U  <J. 

X. 

rg  rsj  2 

•    •>  <  JC 

«-•  »H  CJ 

r-    X   LU  U. 
ZZ~X.CC 

1-1    <L   r-  LU 

Q£  Hm  > 

a.  <  cu  o 


LU 

>  r- 

o    •> 
z 

QC  X 

«l  z 

CU  < 

-J  H- 

CJ  <X 

a  x 

-J   H- 

^  o 
o  u 
u  <i 

a  lu 
cu  a 

co  »— 

<x  •> 

LU   >— < 

cc  a 

CJ 

Z  LU 

»-H    CO 

>  >• 

h-  < 

•-«  a. 

CJ 

a.  u 

<i  —> 

CJ  H- 
CJ 

z  z 

LU  -J 

r-  U- 

ll.  a 

CJ  LU 

z 

cc  <— 

LU  CC. 

-t  LU 
CJ  CJ 

0.  z 

z> 

CU 

X  z 


• 
Q 
LU 

a 
a, 
1-4 

CO 

LU 
CO 


Jt  — ' 


CO 

o 


X 

z 

<l 


X 

X 

o 

CJ 


o 

t— 

CO 

LU 

ac 

LU 

X 


CJ 

•  cc 

CU  CJ 


CJ 

cc 

<l 

CJ 


CJ  i-4  1-* 

a 

-j  z  z 

U   X  (\l  N  U  X  CM 

cc  •»    •»  • 

CJ  CJ  «-•  •-•  CJ  fH 

z  z 

»  1-.  r-   X  "-"  r- 

LU   Q  Z  r-  CJ  Z 

CJ  <.  t-  O  <  «-" 

O    LU  CC   CJ   LU  CC 

-J  X  Q.  <  X  Q. 


* 
* 


* 
* 
* 


* 
* 


too 


3 

10 


U 


I 


I 


CM 

B 

•H 

LX, 


-23- 


PKECISE 

REMARK 

REMARK 
REMARK 
REMARK 
REMARK 
TITLE 
SOLVE 
D1.0*0. 
D2.0833 
DO. 5*0. 
D1.2833 
DO. 3333 
DO. 95 
DO. 25,0 
DO. 7595 
PRECISE 
SOLVE  T 
D1.0*0. 
D0.5»0. 
DO. 3333 
PRECISE 
TITLE 
SOLVE  T 
Dl.0,0. 
D2.45 
DO.  5*0. 
DO. 1428 
DO. 3333 
DO. 1428 
DO.  250, 
DO. 1250 
DO. 20,0 
DO. 1111 
DO. 1666 
DO. 1111 
DO. 1337 
STOP 


THIS  IS  PART  OF  A  TEST  OF  THE  EQUATION 

SOLVER  PORTION  OF  PRECISE  TO  CHECK  THE 

INFLUENCE  OF  THE  MULTIPLE  PRECISION 

OPERATIONS  WHEN  APPLIED  TO  CERTAIN 

ILL-CONDITIONED  MATRICES 

TEST  OF  OPERATIONS  ON  THE  HILBERT  MATRIX  1=4 
THE  FOLLOWING  SYSTEM  OF^EQUATIONS 
5,0.333333333  333333  3  33  33,0.2  5 
333333333333333 
33333333333333333333,0.25  0.2 
33333333*33^3333 
3333  333333333  333,0.25,0.2,0.16666666666666666667 

.2,0.16666666666666666667,0.1428571428571428  5714 
238095238095238095 

TEST  OF  OPERATIONS  ON  THE  HILBERT  MATRIX  1=3 
HE  FOLLOWING  SYSTEM  OF  ©EQUATIONS 
5,0.3333333  33  3  33333  333  33,1.833333  3333333333333 
3333  3  333333  33  333333  3,0.25,1.083  3333333333333333 
3333  3333333  33333,0.2  5,0.2,0.783333333333  33333333 


TEST  OF  OPERATIONS 
HE  FOLLOWING  SYSTEM 
5,0.-333333333333333 

3333333333333333333 
57142857142857143*1 
33333333333333333,0 
5714285714285714,0. 
0.2,0.1666666666666 
»0. lllliillllllllll 
.166  666666666666666 
11111111111111111,0 
6666  666666666  6667,0 
1111111111111111,0. 
662337662^376623 


d)H 


ON 
OF 
33333, 


3,0.25 
.59285 
.25,0. 
125,1. 
666666 
1111*0 
67,0.1 
.1,0.0 
.14285 
1,0.09 


E  HILBERT  MATRIX  1=6 

LINEAR  EQUATIONS 

0.25,0.2,0.166666666666666666667 


,0.2*0.16666666666666666667 

7142857142857143 

2,0.16666666666666666667 

21785714285714285714 

7,0.1428  571428  5714285714 

.99563692063692063692 

428571428571428  5714,0.1251 

84  563692063692063692 

714285714285714*0.125 

09090909090909090909 


Fig.  3.  Typical  input  for  solving  a  set  of  linear  equations. 

The  D  in  the  first  column  of  the  data  cards  is  required  for 
double  precision  input. 
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Table  3 


9° 

tan  0 

cot  e 

1 

0.01745 

50649 

28217 

58576 

51289  0 

57.28996 

16307 

59424 

68727 

815 

89 

2 

0.03492 

07694 

91747 

73050 

04026  3 

28.63625 

32829 

15603 

55075 

651 

88 

3 

0.05240 

77792 

83041 

20403 

88058  2 

19.08113 

66877 

28211 

06340 

675 

87 

4 

0.06992 

68119 

43510 

41366 

69210  6 

14.30066 

62567 

11927 

91012 

805 

86 

5 

0.08748 

86635 

25924 

00522 

20186  7 

11.43005 

23027 

61343 

06721 

086 

85 

6 

0.10510 

42352 

65676 

46251 

15024 

9.51436 

44542 

22584 

92968 

3971 

84 

7 

0.12278 

45609 

02904 

59113 

42311 

8.14434 

64279 

74594 

02382 

5661 

83 

8 

0.14054 

08347 

02391 

44683 

81177 

7.11536 

97223 

84208 

74823 

0566 

82 

9 

0.15838 

44403 

24536 

29383 

88831 

6.31375 

15146 

75043 

09897 

9464 

81 

10 

0.17632 

69807 

08464 

97347 

10904 

5.67128 

18196 

17709 

53099 

4418 

80 

11 

0.19438 

03091 

37718 

48424 

31942 

5.14455 

40159 

70310 

13472 

3221 

79 

12 

0.21255 

65616 

70022 

12525 

95917 

4.70463 

01094 

78454 

23358 

6235 

78 

13 

0.23086 

81911 

25563 

11174 

81456 

4.33147 

58742 

84155 

54554 

6168 

77 

14 

0.24932 

80028 

43180 

69162 

40399 

4.01078 

09335 

35844 

71634 

5715 

76 

15 

0.26794 

91924 

31122 

70647 

25537 

3.73205 

08075 

68877 

29352 

7446 

75 

16 

0.28674 

53857 

58807 

94004 

27581 

3.48741 

44438 

40908 

6506' 1 

6224 

74 

17 

0.30573 

06814 

58660 

35573 

45420 

3.27085 

26184 

84140 

86530 

8856 

73 

18 

0.32491 

96962 

32906 

32615 

58714 

3.07768 

35371 

75253 

40257 

0291 

72 

19 

0.34432 

76132 

89665 

24195 

72658 

2.90421 

08776 

75822 

80257 

9326 

71 

20 

0.36397 

02342 

66202 

36135 

10479 

2.74747 

74194 

54622 

27876 

1664 

70 

21 

0.38386 

40350 

35415 

79597 

14484 

2.60508 

90646 

93801 

53625 

8412 

69 

22 

0.40402 

62258 

35156 

81132 

23481 

2.47508 

68534 

16295 

82524 

0013 

68 

23 

0.42447 

48162 

09604 

74202 

35321 

2.35585 

23658 

23752 

83393 

9587 

67 

24 

0.44522 

86853 

08536 

16392 

23670 

2.24603 

67739 

04216 

05416 

3321 

66 

25 

0.46630 

76581 

54998 

59283 

00062 

2.14450 

69205 

09558 

61635 

6261 

65 

26 

0.48773 

25885 

65861 

42277 

31111 

2.05030 

38415 

79296 

21689 

9011 

64 

27 

0.50952 

54494 

94428 

81051 

37069 

1.96261 

05055 

05150 

58230 

4640 

63 

28 

0.53170 

94316 

61478 

74807 

59159 

1.88072 

64653 

46332 

01236 

0838 

62 

29 

0.55430 

90514 

52768 

91782 

07631 

1.80404 

77552 

71423 

93738 
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0.57735 
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89625 
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91488 
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68877 
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7446 

60 

31 

0.60086 

06190 

27560 

41487 

86644 

1.66427 

94823 

50517 

91103 

0496 

59 

32 

0.62486 

93519 

09327 

50978 

05108 

1.60033 

45290 

41050 

35532 

6733 

58 

33 

0.64940 

75931 

97510 

57698 

20629 

1.53986 

49638 

14582 

90482 

6797 

57 

34 

0.67450 

85168 

42426 

63214 

24609 

1.48256 

09685 

12740 
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7157 

56 

35 

0.70020 

75382 

09709 

77945 

85227 

1.42814 

80067 

42114 
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0618 

55 

36 

0.72654 

25280 

05360 

88589 

54668 

1.37638 

19204 

71173 

53820 

7210 

54 

37 

0.75355 

40501 

02794 

15707 

39564 

1.32704 

48216 

20410 

03715 

9473 

53 

38 

0.78128 

56265 

06717 

39706 

29500 

1.27994 

16321 

93078 

78031 

1030 

52 

39 

0.80978 

40331 

95007 

14803 

69914 

1.23489 

71565 

35051 

39855 

6175 

51 

40 

0.83909 

96311 

77280 

01176 

31273 

1.19175 

35925 

94209 

95870 

5308 

50 

41 

0.86928 

67378 

16226 

66220 

00956 

1.15036 

84072 

21009 

55587 

6331 

49 

42 

0.90040 

40442 

97839 

94512 

04772 

1.11061 

25148 

29192 

87014 

3482 

48 

43 

0.93251 

50861 

37661 

70561 

21856 

1.07236 

87100 

24682 

53294 

6028 

47 

44 

0.96568 

87748 

07074 

04595 

80273 

1.03553 

03137 

90569 

50695 

8833 

46 

45 

1.00000 

00000 

00000 

00000 

00000 

1.00000 

00000 

00000 

00000 

0000 

45 

cot  0 

tan  0 

0° 

Values  of  the  two  coefficients  <r4  and  a6  and  the  three  function*  « 

For  the  convenience  of  the  reader,  there  is  included  in  Table  3  a  compilation 
of  25D  values  of  the  tangent  and  cotangent  for  arguments  1°(1°)89°.  These  data, 
which  are  required  in  computing  the  values  of  c,  are  here  tabulated  with  the  same 
range  and  precision  as  for  the  values  of  sine  and  cosine  given  by  G.  W.  and  R.  M. 
Spenceley  [4].  The  only  comparable  table  of  decimal  approximations  to  the  tangent 
appears  to  be  the  relatively  inaccessible  30D  table  of  Herrmann  [5], 


Fig.  6.  See  the  following  figures  for  the  degree  to  which  this  table 
layout  is  approximated  via  the  instructions  on  the  opposite 
page. 
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Figure  7.  Results  from  the  program  shown  in  Figure  5. 
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PRECISE 

TITLE  ELEMENTARY  AND  TRANSCENDENTAL  FUNCTIONS  TO  25  FIGURES 

FORMAT  ON  1  CARC 

1C  1F2.9,1P5K1S7(1E25.45  » 

GENERATE  -1.{.C5)1.  IN  CCL  1 

ARCSIN  1,2 

HEADING  ON  1  CARC 

X  ARCSIN(X) 

PRINT  1,.2 
ARCCCS  1,3 
HEADING  CN  1  CARC 

X  ARCCGS(X) 

PRINT  1,3 

GENERATE  iCK  .CI)  .  1  (  .1 )  1.  (1.  )  1C.  ( 10.  )  100.  ( 100.  UOOO.  IN  COL  4 

ARCTAN  4,5 
HEADING  CN  1  CARC 

X 
PRINT  4,5 
ARCCCT  4,6 
HEADING  CN  1  CARC 

X 
PRINT  4,6 
FIXED  WITH  28 
HEADING  CN  1  CARC 

ARCSIN(X) 
PRINT  2 
HEADING  CN  1  CARC 

ARCCCSU) 
PRINT  3 
HEADING  CN  1  CARC 

ARCTAN(X) 
PRINT  5 
HEADING  CN  1  C/SRC 

ARCCCT(X) 
PRINT  6 

FLCATING  WITH  28 
HEADING  CN  1  CARC 

ARCSIN(X) 
PRINT  2 
HEADING  CN  1  CARC 

ARCCCS(X) 
PRINT  3 


ARCTAN(X) 


ARCCCT(X) 


Figure  8.  A  portion  of  a  PRECISE  program  to  compute  tables  of 

trigonometric  functions.  See  Tables  8,  9,  and  10  for  some 
results  of  this  calculation. 
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Figure  11.  A  portion  of  a  table  of  trigonometric  function 
computed  on  PRECISE.  The  values  v/ere  checked  against 
existing  tables  up  to  the  snaded  area. 


-32- 


at 


Q. 


LJ 


OO 


3t 

U> 

l_) 

at 
u. 

tu 
o 
at 

_; 

LJ 


00 

I— 

cj 

LU 
Ot 


o 
lj 


Q£ 


o 


< 

UL 
UJ 

UJ   (_J 


in 

* 


01 


o   IN 


lj»  at 
to  a. 


r\l  00 


-C  _3 


(N 


*i  «a. 


UL 


a 

at 
<i 
ji 
at 
a 


LJ 

at 
a 
<_> 

UJ 

at 

CM 


XJ 

Ul 

2. 

Z5 


CM 

UJ 

J. 


•<=.    A. 


>-«  at 


oo 


at 

<». 

JE 
UL 
U 


I"- 
"O 

to 

00 


LJ   UJ 


0O 

LJ 
LL 


UJ 


LJ    0O 

h-   £- 
>    UJ 


:>  t— 

UJ   LO 
UJ    LU 


<£   LL 

oo 

<£. 
<l 
at 


m 


UJ 


<i 

UL 


UJ 


en 


3 

Ul. 


00 


LU 

at 

X 


-~  LU 
LU  U- 
UJ   < 

ot  *- 

x  — 

•—  a 

x 

LU   <t 

a.  -k 
<:  at 

<~   LJ 


LJ 

UJ 


LJ 

at 

<i 

JE 

it 

aj 

LJ 

at 


•■  at 


cu 

sU    LJ 


V>J 


LJ 
O 

LU 
OO 

«1 

at 

LU 


LU  LU 

LU  LU 

ot  at 

x  x 


LU  LU 

a.  a. 


Jt 


(M 


-i. 

LJ 

at 

LU 


at 


UJ  LU 

LU  LU 

at  ot 

x  x 


LU  LU 

a.  cw 
<t  <t 


CO 


t\J 


>jn  t— 

O 

I—  «c 

<c  <— 
I— (  Jfc 
at  lu 
a.  at 


LU 

LU 

at 
x 


LU 

a. 


u. 


at 

<i 
-j 

LJ 


LJ 


aj 


-J  _J 

LU    —I 
^t    LJ 

—I   LJ 

X 

LJ   I— 

r—   <. 

— •     1-4 

jc  at 
00  u. 


00 


W 
CO 
H 

« 

c 

•H 

W 

aj 
u 

+j 
to 

cd 
iw 

d 
o 

-H 
-U 
rd 

rH 

Q< 
-H 

a 


(D 

a 

(0 

-P 

Q) 
-C 
-M 

MJ. 
0 

CO 

CD 

rH 

s 
m 
x 
a) 

0) 

e 
o 
co 


CN 


CD 

■H 

fa 


33- 


APPENDIX  I  -  JOB  DEFINITION  OF  PRECISE  PROGRAMS 


PRECISE  was  written  to  operate  as  a  sub-monitor  under  the  IBSYS  or 
DC-IBSYS  Monitor  System  for  the  IBM  7090/7094,  or  as  an  IBMAP  program  to 
be  run  under  the  IBJOB  monitor  under  IBSYS.     Either  of  the  two  versions 
may  be  produced  by  setting  one  parameter  and  making  an  absolute  or  re- 
locatable assembly  of  the  source  program  via  the  IBMAP  Assembler  of  IBJOB. 

JOB  Definition  when  PRECISE  is  a  Monitor 


A  PRECISE  job  will  consist  at  most  of  the  following  cards: 

$JOB 

$EXECUTE  PRECISE 

$ID 

[  SOURCE  1  ] 

[  SOURCE  n  ] 
STOP 
end  of  file 

Where  [  SOURCE  i  ]  consists  of  the  PRECISE  command  followed  by  as  many 
as  desired  of  the  commands  in  the  PRECISE  language  (excluding  the  PRECISE 
and  STOP  commands). 

The  requirements  as  to  need  and  format  of  the  $JOB  and  $ID  cards  will 
vary  among  installations. 

JOB  Definition  with  PRECISE  under  IBJOB 


A  PRECISE  job  when  using  the  relocatable  deck  under  the  IBJOB 
Monitor  will  be  made  up  as  follows. 


$JOB 

$EXECUTEAAAAAAAlBJOB 
$ID 
$ IBJOB 

[  PRECISE  program  decks  ] 
$DATA 

[  SOURCE  1  ] 


[  SOURCE  n  ] 
STOP 
end  of  file 

Where  [  SOURCE  i  ]  consists  of  the  PRECISE  command  followed  by  as 
many  as  desired  of  the  commands  in  the  PRECISE  language  (excluding  the 
PRECISE  and  STOP  commands). 
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APPENDIX  II  -  FORMAT  DESCRIPTION 

A  format  consists  of  any  number  of  tield  specifications. 
The  general  field  specification  is  of  the  following  form. 

aSb  (c) gPnTd«wZ 

a,b,c,g,n,d,  and  w  are  positive  integers. 

T  is  the  type  conversion  character  and  T  is  I,  E,  F,  0, 
A,  or  X. 

Z  is  a  slash,  comma,  or  blank,  and  Z  serves  as  a  field 
specification  separator. 

Each  non-numeric  Hollerith  character  in  a  field  specifi- 
cation indicates  what  is  to  be  done  with  the  integer  just 
preceding  it. 

Each  field  specification  is  executed  n  times  and  n  data 
fields  are  printed,  each  having  a  total  column  width  w  with 
the  spread  of  information  being  determined  by  a,b,c,g,  and 
d. 

T  =  I  specifies   full   word   binary   to   decimal   integer 
conversion, 

T  =  E  specifies   binary  to  floating  point  decimal  conver- 
sion of  multiple  precision  numbers. 

T  =  F  specifies  binary  to  fixed  point  decimal   conversion 
of  multiple  precision  numbers. 

T  =  0  specifies  full  word  binary  to  octal  conversion. 

T  =  A  specifies  BCD  to  BCD  conversion. 

T  =  X  specifies  insertion  of  blank  characters. 

If  either  of  the  integers  a,b,c,  or  g  is  zero  then  that 
integer  and  the  following  Hollerith  character  may  be  omit- 
ted. 

The  format  is  scanned  from  left  to  right.  Couversion  as 
indicated  by  a  specification  is  completed  before  checking  on 
the  separator  Z.  If  Z  is  a  blank,  then  all  information 
converted  at  this  point  is  printed  and  if  there  is  still 
more  information  to  be  converted  and  printed,  scanning  goes 
back  to  the  beginning  of  the  format.  If  Z  is  a  slash,  all 
information  converted  at  this  point  is  printed  and   scanning 
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for  the  next  line  starts  immediately  alter  the  slash.  K 
consecutive  slashes  separating  two  specifications  cause  K- > 
blank  lines  to  be  printed.  If  1  is  a  comma,  conversion 
continues  with  the  specification  following  the  comma  being 
scanned. 


Binary  to  _  Decimal  Integer  ^Conversion [T=  1} l3.RH.LlLl. 

n  single  celled  binary  integers  are  converted  and  printed 
as  decimal  integers,  each  with  a  total  column  width  w;  and 
each  number  is  enclosed  in  parenthesis  if  g#0.  The  restric- 
tions are  n>0,  w<26. 

Example:  Suppose  the  binary  equivalent  of  the  integers  1, 
-2,  3,  7000,  and  -56789  are  to  be  printed.  Then  the  format 
1P1I4,4I5  would  cause  the  following  line  to  be  printed. 

(1)    -2     3  700056789 

Note  that  the  negative  sign  of  the  last  number  would  not  be 
printed  since  w  =    5. 

Bin ary_to_ Floating  Point  Decimal  .Conversion 
JT=  E] JaS  bi'cjjjPnEd  •  w  zj 

n  3  cell  multiple  precision  numbers  are  converted  and 
printed,  each  number  having  a  total  column  width  w.  Infor- 
mation within  the  column  is  arranged  as  follows  from  left  to 

right. 

1.  w-a-b-c-d-2  blanks. 

2.  the   decimal  exponent   of  the  number  occupying  b  places 
and  in  parenthesis  if  g*0. 

3.  The  sign  of  the  number  (blank  for  *,  or  a  minus  sign) . 

4.  a  significant  digits  of  the  number. 

5.  a  decimal  point. 

6.  d  significant  decimal  digits  (rounded). 

7.  the  decimal  exponent  of  the  number   occupying   c   places 
and  in  parenthesis  if  g*0, 

Restrictions  are  n>0, w>a+b+c+d+2 ,  b<26,  c<26,  d+a<30. 
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Example:.  The  format  1S4 ( 1P  1E5.  13, 3)  1E10. 19  would  cause  the 
binary  equivalent  of  12345.67  and  -.00765432  to  be  printed 
as  follows. 

(4)  1.23457     -.7654320000  -2 

Binary  to  Fixed  Point  Decimal  ConversionjT^Fj faSnFd«wZ] 

n  3  cell  multiple  precision  numbers  are  converted  and 
printed  as  fixed  point  decimal  numbers,  each  number  having  a 
total  column  width  w.  Information  within  the  column  is 
arranged  from  left  to  right  as  follows.  Before  conversion, 
the  number  is  multiplied  by  10  to  the  a  power. 

1.  The  integral  part  of  the  number   (with   sign)   occupying 
w-d-1  spaces. 

2.  A  decimal  point. 

3.  d  decimal  places  (rounded)  of  the  fractional  part  of  the 
number. 

fiestrictions   are   n>0,   w>d+2+number   of  digits  in  integral 
part. 

This  subroutine  arbitrarily  refuses  to  print  fixed  point 
multiple  precision  numbers  whose  absolute  value  is  greater 
than  23S-1.  If  such  a  number  is  found,  the  specification  is 
treated  as  E  type  rather  than  F  type  and  the  number  is 
printed  in  floating  point  form,  but  the  power  of  ten 
exponent  cannot  be  printed  unless  either  b (or  c)  is  written 
in  the  specification.  b  and  c  are  ignored  in  F  type 
conversion,  so  no  harm  is  done  by  specifying  one  of  the  two 
just  in  case  large  numbers  are  printed. 


Example:  The  format  3J3F10.20  would  cause  the  binary 
equivalent  of  -101.01,  10*o,  and  .00004778966  to  be  printed 
as  follows. 

-101.0100000000       .1000000000  41  .0000477897 

Binary  to  Octal  Conversion  (T=0)   faSb (c) gPnOd* vZ] 

n*  single  celled  binary  integers  are  converted  and 
printed  as  octal  integers,  with  locations  if  desired.  Each 
number  has  a  total  column  width  w  and  information  within  the 
column  is  arranged  from  left  to  right  as  follows. 
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1.  w-a-b-c-d  blanks. 

2.  The  location  as  a  decimal  integer  occupying  b  places  and 
in  parenthesis  if  g  =  1,  3,  9,  or  11. 

3.  The   location  as  an  octal  integer  occupying  c  places  and 
in  parenthesis  if  g  =  2,  3,  10,  or  11. 

4.  a  blanks. 

5.  The  last  d  octal  digits  of  the  number. 

Restrictions  are  d<12,  b<26,  c<26,  v>a+b+c+d,  n>0.  If  g>8, 
then  as  far  as  the  calling  seguence  is  concerned,  the  number 
described  above  is  not  counted  as  having  been  printed;  thus 
allowing  a  single  cell  to  be  printed  both  in  octal  and 
decimal  with  or  without  decimal  and/or  octal  locations;  or 
if  d  =  0  allowing  a  location  to  be  printed  with  either  a 
single  celled  number  or  multiple  precision  number. 

Example:  The  format  1S5  (9P1oO. 7, 116, 2S3)  1o12. 21 , 1o3. 7  would 
cause  the  octal  numbers  -17,314631463146  and  7777777  in 
decimal  locations  100,  101,  and  102  to  be  printed  as 
follows. 

(100)    -15     145   314631463146     777 

*  If  n>1  then  g  must  be  less  than  8. 

BCD  Conversion  {T=A) fnAwZ] 

n  BCD  words  are  printed  as  n  BCD  fields  with  each  field 
having  a  total  column  width  w.  If  w>6 ,  then  each  field 
consists  of  w-6  blanks  followed  by  the  BCD  word  consisting 
of  6  characters.  If  w<6  then  the  field  consists  of  the 
leftmost  w  characters  of  the  BCD  word. 

Blank  Insertion  JT=X) InXH 

n  blanks  are  inserted  into  the  print  line. 

BCD  C ha racter§_ Within  .the  Format 

All  of  the  previously  described  field  specifications  may 
be  preceded  by  iC  followed  by  i  Hollerith  characters.  The  i 
characters  will  be  printed  just  before  the  first  number 
printed  by  the  specification  which  follows. 
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Printing  of  Tables fmKaSbfc) qPnTd«zW} 

To  make  numbers  with  many  digits  more  readable,  an 
additional  feature  is  provided  in  the  format  specification 
for  E  and  F  type  conversion.  Either  the  F  or  the  E 
specification  as  described  previously  may  be  preceded  by  mK. 

The  effect  of  mK  is  to  cause  blocks  of  ra  digits  to  the 
right  of  the  decimal  point  to  be  separated  by  one  blank. 
The  separation  will  occur  on  both  E  and  F  type,  and  it  will 
continue  until  reset  by  another  mK  appearing  in  some 
specification. 

When  this  feature  is  used,  the  column  width  w  must  be 
increased  enough  to  allow  j  extra  spaces. 


Let  d/m=k*fi,  k  an  integer  and  0<B<m. 
Then  j  = 


{0  if 
k  if 
k-1  i 


integer 

m=0. 

R>0. 

f    R=0. 
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APPENDIX  III  -  SYSTEMS  INFORMATION 

The  PRECISE  system  consists  of  about  8000  IBMAP  instructions  in 
the  form  of  a  PREST  deck,  which  is  about  3/4  of  a  box  of  cards.  To 
assemble  PRECISE  as  an  absolute  program  (which  can  be  edited  into  IBSYS) 
use  the  ALTER  feature  of  IBJOB  to  replace  card  2  with  the  IBMAP  symbolic 
instruction: 

RELMODASETAAAAA2 

and  use  the  ABSMOD  option  on  the  $  IBMAP  card.  To  make  PRECISE  a  sub- 
monitor  under  IBSYS,  enter  "PRECIS"  into  the  name  table  of  IBSYS,  and 
then  insert  the  absolute  program  deck  into  the  desired  position  of  the 
system  tape. 

A  small  dummy  IBMAP  program  is  included  with  the  system  and  this 
program  should  be  discarded  if  making  an  absolute  assembly. 

To  assemble  PRECISE  as  a  relocatable  program,  use  the  ALTER  feature 
of  IBJOB  to  replace  card  2  with  the  IBMAP  symbolic  instruction: 

RELMODASETAAAAAl 

and  include  the  small  dummy  IBMAP  program  to  be  assembled.  This  dummy 
program  has  the  deck  name  USERS  and  has  the  entries  USERSA,  USERSB,  USERSC, 
USERSD,  and  USERSE.  This  program  is  necessary  because  a  relocatable  assembly 
of  the  system  generates  calls  to  the  above  entries.  These  external  calls 
are  automatically  deleted  from  an  absolute  assembly.  The  USERS  program 
also  has  a  common  block  definition  of  key  calls  of  the  system.  A  complete 
listing  of  USERS  is  given  later.  The  relocatable  assembly  of  PRECISE  also 
defines  many  entry  points  to  various  parts  of  the  system.  By  use  of  the 
entries  and  COMMON  values,  new  commands  may  be  implemented  in  the  system 
when  the  code  to  accomplish  the  command  function  is  written  and  made  part 
of  USERS  program.  Ways  of  adding  commands  are  given  later. 

The  PRECISE  system  consists  of  5  blocks  of  coding. 

1.  The  system  control  part  of  the  system. 

2.  UOMAMPP  package,  modified  for  IBMAP  and  I/O  calls  to  UOMAIOS  (3). 

3.  UOMAIOS  package  by  A.  Beam  and  available  as  SHARE  distribution  number  3444. 

4.  The  PRECISE  coding. 

5.  The  vocabulary  of  commands. 

Available  storage  is  computed  by  using  .JCOR  in  relocatable  versions 
and  SYSCOR  in  absolute  versions. 

The  entry  points  in  relocatable  versions  are 

1.  All  entries  defined  in  the  UOMAMPP  write-up. 

2.  All  entries  defined  in  the  UOMAIOS  write-up. 

3.  INVLD  which  may  be  called  when  an  invalid  command  is  detected  by  the 
user's  program. 
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The  Vocabulary 

Each  legal  command  defined  to  the  system  must  have  a  two  word  entry 
in  the  vocabulary.  Only  two  types  of  entries  are  described  here,  and 
should  be  sufficient  for  expansion  of  the  system.  The  entry 

BClAAAAl,  name 
TXLAAAAcode,4,0 

is  used  when  the  coding  at  location  code  is  part  of  the  main  system. 

The  entry 

BClAAAAl,  name 

CALLAAAcode 

0RGAAAA*-2 

is  used  when  the  coding  at  location  code  is  external  to  the  system.  Name 
is  the  PRECISE  command  word  and  consists  of  6  or  less  alphabetic  characters, 
right  justified  with  leading  zeros.  For  example,  the  command  word  SQRT 
appears  in  the  first  entry  type  above  as  BCI   1,00SQRTA.  The  coding  at 
location  code  must  be  a  subroutine  which  is  entered  via  a  TSX   code,  4  and 
which  exits  via  a  TRAAl,  4A.  The  coding  utilizes  the  COMMON  block,  and 
entry  points  of  PRECISE  to  carry  out  the  desired  function.  Entries  to  the 
vocabulary  must  be  in  ascending  sort  order  on  the  BCI  command  name  words. 

The  COMMON  block 

A  listing  of  the  COMMON  block  appears  in  the  listing  of  USERS  program. 
To  explain  the  most  used  values  in  the  block,  suppose  the  following  PRECISE 
command  has  been  received  by  the  system. 

name   C(l),  C(2),...,  C(n) 

Where  name  is  the  command  word  and  C(i),  i=l,...,n  are  arguments  to  the 
command. 

When  control  is  given  to  the  code  for  execution  of  the  command,  the 
following  information  is  available  in  the  COMMON  block. 

S    :  The  address  portion  of  S  contains  the  matrix  origin. 

NC   :  The  address  portion  of  NC  contains  the  maximum  column  number  available. 

NR   :  The  address  portion  of  NR  contains  the  maximum  row  number  available. 
Hence  the  current  size  of  the  matrix  is  the  product  of  NR  by  NC. 

ARGS  :  The  origin  of  a  vector  containing  the  arguments  C(i)  of  the  command, 

3  words  per  argument.  Argument  C(i)  will  be  in  locations  ARGS+3(i-l)+j, 
j»0,l,2;  if  C(i)  is  a  multiple  precision  number;  and  in  location 
ARGS+3i-l  if  C(i)  is  an  integer. 
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STATS       :     The  origin  of  a  vector  containing  one  word  per  C  (i) .     The 
word  at  STATS+i-1  will  be  0  if  argument  C(i)  is  an  integer, 
and  will  be  non-zero  if  argument  C(i)  is  a  multiple  precision 
number. 

PICKUP     :     The  origin  of  a  vector  containing  one  word  per  C(i).     Each 
word  at  PICKUP+i-1  is  of  the  form 


PZE      x,I,A 

where  x  may  be  any  value. 

I  will  be  0  if  C(i)  is  a  multiple  precision  number. 
I  will  be  1  if  C(i)  is  an  integer. 

A  is  the  origin  of  the  vector  C(i)  if  1=1  or  the  origin 
of  a  multiple  precision  number  if  1=0. 

NARGS   :  The  address  portion  of  NARGS  contains  n=  the  number  of  arguments 
on  the  PRECISE  command  currently  under  consideration. 

NRA    :  The  address  portion  of  NRA  is  the  current  operating  length  of 
columns  of  the  matrix.  Operations  on  rows  should  be  carried 
out  on  the  1st  through  NRA  th  row. 

BUFFER  :  The  PRECISE  command  (input  card  image)  currently  under  consideration 
is  stored  in  RJFFER+I,  1=1,2,.. . ,12. 

EXTRAS  :  The  origin  of  a  100  word  vector  which  may  be  used  as  temporary 
storage. 

The  above  information  is  sufficient  for  a  subroutine  to  operate  on  one  or 
more  columns  in  any  desired  fashion.  Columns  should  be  operated  upon  only 
through  the  NRA  th  row  of  the  matrix. 

Expanding  the  System 

It  may  be  desirable  to  add  commands  to  the  PRECISE  language.  Coding  for 
new  commands  should  probably  be  checked  out  using  the  relocatable  IBJOB  version 
of  PRECISE. 

The  easiest  way  to  add  and  check  out  a  new  command  is  to  use  one  of  the 
entries  in  the  USERS  program.  Then  at  a  later  time  the  permanent  vocabulary 
entry  and  the  coding  can  be  made  a  part  of  the  system. 

A  listing  of  the  USERS  program  follows,  and  includes  the  COMMON  block 
which  is  necessary  for  any  new  relocatable  programs  which  may  be  added  to  the 
system. 
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SIBMAP  USERS    ( )0K 


*  USERS  PROGRAM  TO  ADD  AND  CHECK  FUNCTIONS 
ENTRY    USERSA 


ENTRY    USERSB 
ENTRY    USERSC 


ENTRY    USERSD 
ENTRY    USERSE 


*  C  0  M  M  0  N    BLOCK 


* 


USE      // 


S       BSS      1   ADDRESS  DEFINES  MATRIX  ORIGIN 

NC_     BSS 1   MAXIMUM  NUMBER  OF  COLUMNS  AVAILABLE 

"NMAX    BSS      1   NUMBER  OF  CELLS  IN  MATRIX 

NR BSS 1   MAXIMUM  NUMBER  OF  ROWS  IN  MATRIX 

"N"R"7\  B~S~S"  I       COUNT    OF    ROWS    READ     IN 

SKIPOP  BSS 1   0  UNTIL  AN  ID  IS  FOUND 

EXTRAS  BSS      100  USED  TO  STORE  SUMS AVAILABLE  FOR  TEMP.  USE 

SUMSWT  BSS 1   NOT  0  IF  COLUMN  SUMS  WANTED 

TYPECH  BSS      1   TYPE  OF  PR  I NT  =  0 , 1 , . . . » 10 

BLOCKC  BSS 1   PRINT  BLOCK  COUNT 10  LINE/BLOCK 5  BLOCKS/PAGE 

ERSKIP  BSS      I   SET  NOT  0  WHEN  BAD  OP  ENCOUNTERED 

I NVEC   B  SS 1   SET  NOT  0  BY  'READ*  COMMAND 

NERROR  BSS      1   ERROR  TYPE  IN  DECREMENT  =  0»1»2»... 

LINENO  BSS 1   LINE  COUNT  FOR  PRINTING 

PAGENO  BSS      1   PAGE  COUNT  FOR  PRINTING 

PICKUP  BSS 36  PZE  XtltA 1=0  IF  CONSTANT , ELSE  1 A  =  ARG  ADDRESS 

ARGS    BSS      3*36  START  OF  ARGUMENT  VECTOR INTEGERS  ARE  IN  3RD  CELL 

STATS   BSS 36  STATS  VECTOR INTEGER  IF  QtELSE  MPP  NUMBER 

BUFFER  BSS      21  INPUT  CARD  IMAGE  IN  BUFFER+I .  I=l*2»... 

NARGS   BSS 1   NUMBER  OF  ARGUMENTS  ON  CARD 

ENDMAT  BSS      1  ADDRESS  DEFINES  LAST  CELL  AVAILABLE  FOR  MATRIX 

CONTRL   // 

USE      PREVIOUS 

# 

"*  END  OF  COMMON  BLOCK 

*  USE  SYMBOLS  J»H(LESS  THAN  6)»K(EXCEPT  K ) » Q » V » X » Y ( ALL) 

USERSA  NULL     * 

USERSB  NULL     * 

USLRSC  NULL     * 


USERSD  NULL     * 


USERSE  NULL     * 
TR/\ 1,4 

END 
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Note  that  all  the  entries  in  the  USERS  program  do  nothing  except 
return  to  the  main  PRECISE  package. 

Hence,  one  could  write  the  commands 

USERSA 


USERSE 

with  any  arguments  desired.  The  commands  would  be  executed,  but  nothing 
useful  would  be  accomplished.  However,  on  entry  to  the  USERS  program 
certain  information  would  have  already  been  set  up  in  the  COMON  block; 
namely  the  PICKUP,  ARGS,  and  STATS  vectors  which  define  the  arguments  on 
the  command  and  NARGS  which  gives  the  number  of  arguments.  This  infor- 
mation along  with  NRA  (which  describes  the  current  number  of  rows  in  the 
matrix)  is  sufficient  to  operate  on  the  arguments  which  appear  on  the 
command. 

For  example,  suppose  we  need  a  PRECISE  command  which  computes  the 
product  of  all  the  numbers  in  the  column  specified  by  the  first  argument 
and  stores  this  product  in  the  column  specified  by  the  second  argument 
of  the  command.  For  check  out  we  could  call  the  command  name  USERSA  and 
by  writing  the  command 

USERSA  COL  5,4 

the  product  of  all  values  in  column  5  will  be  stored  over  all  values  in 
column  4.  The  coding  listed  below  will  produce  the  desired  result  if  the 
entire  block  of  coding  is  inserted  into  the  USERS  program  (listed  above) 
between  the  two  cards  labeled  USERSA  and  USERSB,  and  the  USERS  program 
is  re- assembled. 
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*    USERSA    C1.C2    INPUT    CI    AND    OUTPUT    IN    C2 


*  COMPUTE  A(l) 
PRODUC  NULL 

*A(2  )*...*A(NRMAX)  AND  STORE  IN  SPECIFIED  VECTOR 

SXA 
AXT 

UA5,4 
0,1 

LXA 
CLA 

NRA.2 
PICKUP 

STD 
CLA 

UA2 
PICKUP+1 

STD 
TSX 

UA4 
ENTRY, 4 

CLA1,0,UA6 
DPFLTI 

UA1 

TSX 

DPEXIT 
ENTRY, 4 

UA2 

MPY»1»** 
DPEXIT 

TXI 
TIX 

*+l»l»-3 

UA1,2»1 

AXT 
LXA 

0,1 
NRA,2 

UA3 
UA4 

TSX 

ENTRY, A 
ST0>1»** 

TXI 

DPEXI.T 
*+l,l,-3 

UA5 

TIX 
AXT 

UA3,2,1 
**,4 

UA6 

TRA 
PZE 

1,4 

1 

If,  at  some  later, time,  it  is  desired  to  make  the  product  command 
a  permanent  part  of  the  system,  then  a  vocabulary  entry  (for  example, 
PRODUC)  could  be  made  and  the  above  coding  could  be  inserted  into  the 
PRECISE  system.     After  a  re- assembly,  the  command 

PRODUCT  OF  C(l)  TO  C(2) 

would  be  a  part  of  the  PRECISE  language.  The  vocabulary  entry  would 
consist  of  the  two  cards. 

BClAAAAAl, PRODUCT 
TXLAAAAAPRODUC 

which  would  be  inserted  between  the  entries  for  PRNT10  and  RAISEI, 

The  above  method  of  adding  a  command  will  allow  the  command  PRODUCT 
to  work  in  both  the  relocatable  and  absolute  versions  of  the  system. 
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APPENDIX  IV 


MULTIPLE  PRECISION  PACKAGE 

,  (UCM  MPP) 
for  the 
IBM  7090/7094 

by 
Alfred  E.  Beam 


Part  of  the  work  on  this  package  was  done  while  the  author  was  em- 
ployed at  the  National  Bureau  of  Standards.  The  final  version  of  the 
package  and  its  incorporation  under  IBSYS  was  performed  at  the  Computer 
Science  Center  of  the  University  of  Maryland  and  was  supported  by  grant 
NsG-398  of  the  National  Aeronautics  and  Space  Administration. 
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Introduction 

Many  scientific  problems  lead  to  computational  requirements 
involving  extensive  accuracy  or  extrememly  large  data  ranges. 
This  report  describes  an  arithmetic  program  package  for  the  IBM 
7090  or  IBM  7094  computer  which  essentially  extends  the  standard 
word  length  and  the  normal  data  range  for  these  machines .  More 
precisely  the  package  interprets  and  executes  a  set  of  pseudo- 
instructions  in  a  multiple  precision  mode  and  provides  input- 
output  for  several  types  of  data,  in  particular  multiple  precision 
numbers.   Pseudo-instructions  are  provided  not  only  for  the  stan- 
dard arithmetic  operations  but  also  for  most  of  the  elementary 
functions.   In  addition  certain  subroutines  have  been  provided  for 
solving  a  large  class  of  problems  requiring  multiple  precision 
computations . 

Numbers  operated  upon  are  in  a  floating  point  form.   Internally 
a  number  occupies  3  machine  cells;  the  fractional  part  of  the 
number  is  stored  in  the  first  two  cells  and  in  the  first  part  of 
the  third  cell.   The  power-of-two  exponent  is  stored  in  the  re- 
maining part  of  the  third  cell.   The  normal  format  of  the  third 
cell  specifies  that  26  bits  are  for  the  fractional  part  of  the 
number  and  the  last  9  bits  are  for  the  power-of-two  exponent. 
A  subroutine  has  been  included  in  the  package  which  makes  it 
possible  to  reset  the  format  of  the  third  cell  and  thus  vary  the 
possible  range  of  the  data. 

In  brief,  the  package  will  effectively  transform  a  7090 
into  a  machine  which 

1.  is  about  3  1/2  times  as  accurate. 

2.  is  about  20  times  slower. 

3.  has  about  1/3  the  storage  capacity. 

4.  has  the  capability  of  handling  numbers  whose  magnitudes 
are  very  large  or  very  small.   For  example,  the  package  can  compute 
e  to  the  millionth  power  correct  to  22  significant  decimal  digits, 

This  package  has  already  found  wide  acceptance  and  has  been 
used  extensively  at  the  University  of  Maryland,  the  NASA-Goddard 
Space  Flight  Center,  the  National  Bureau  of  Standards,  Washington, 
D.C.,  and  the  General  Motors  Corporation,  Detroit,  Michigan.   Its 
main  uses  have  been  for  computing  tables,  checking  of  double 
precision  IBM  7090  codes  and  for  calculation  of  rational 
approximation  coefficients. 
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Identification 


a.  7090/7094  Multiple  Precision  Package   UOM  MPP 

b.  A.  Beam,  September  23,  1963 

c.  Computer  Science  Center,  University  of  Maryland,  College  Park, 
Maryland 


Purpose 

To  interpret  and  execute  certain  pseudo-instructions  in  a 
multiple  precision  mode,  and  to-  provide  input-output  for  several 
types  of  data  including  multiple  precision  numbers.   Also  to  pro- 
vide a  few  subroutines  so  that  the  package  can  be  used  for  solving 
a  large  class  of  problems  which  require  high  precision  and  or  com- 
putation using  numbers  of  very  great  or  very  small  magnitudes. 


Restrictions 

The  symbolic  deck  for  MPP  contains  many  symbols  and  is  set  up 
as  a  relocatable  FAP  subroutine  with  many  entries .   In  this  form  it 
is  usable  as  a  subprogram  by  programs  which  run  under  the  FORTRAN  II 
Monitor  under  IBSYS.   All  I/O  is  self-contained  except  for  the  use 
of  IOEX,  hence  the  package  could  be  used  under  any  system  which 
operates  under  IBSYS . 

Underflow-overflow  triggers  and  indicators  are  not  saved  by 
the  package. 


Method 

In  the  interpretive  mode  successive  pseudo-instructions  are 
brought  into  the  accumulator  and  control  is  sent  to  the  proper 
place  within  the  package  for  execution.   Numbers  operated  on,  and 
results  given,  while  in  the  interpretive  mode  are  of  the  form 

N=F'2P,  1/2  <  Jf|  <  1  or  F=0,  and  -2b_1  <   P  <   2b_1,35>b2  9 


N  is  stored  in  3  consecutive  cells  with  the  sign  of  F  in  each 
cell.  There  are  b  (at  present  b=9,  but  this  may  be  changed  by  use 
of  a  subroutine  to  be  described  later)  bits  at  the  end  of  the  third 
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cell  which  contain  P+2    and  the  105-b  bits  (not  including  signs) 
of  the  3  cells  contain  JFJ.   Hence  "multiple  precision  number"  in 
this  write-up  is  taken  to  mean  a  number  which  has  at  least  21  but 
no  more  than  28.8  significant  decimal  digits,  and  the  potential 
range  of  a  number  is  increased  as  the  number  of  significant  digits 
is  decreased. 

Arithmetic  operations  are  coded  using  fixed  point  machine 
instructions  and  all  yield  a  (105-b)  bit  rounded  result.   On  under- 
flow the  result  is  automatically  set  to  zero.   If  overflow  is  detected, 
a  message  is  printed  and  execution  is  normally  terminated.   A  subroutine 
is  described  later  which  allows  the  user  to  get  control  if  overflow 
is  detected. 

The  method  used  in  the  elementary  functions  is  described  later. 

Usage 

Usage  of  subroutines  is  described  later.   The  interpretive  mode 
is  entered  by  giving  the  instruction 

L       TSX      ENTRY , 4 

and  location  L  is  followed  by  as  many  as  desired  of  the  pseudo- 
instructions.   Interpretation  begins  at  location  L+l  and  continues 
until  one  of  several  pseudo-instructions  is  found  which  cause  de- 
parture from  the  interpretive  mode. 

Storage 

MPP  requires  about  (4000)    storage  locations.   COMMON  storage 
is  not  used. 


Timing 

Timing  for  a  code  using  MPP  should  be  from  20  to  25  times  that 
for  the  same  code  using  single  precision  floating  point  arithmetic. 

All  of  the  elementary  functions  (except  the  square  root  and 
integer  exponentiation)  in  the  package  are  computed  by  series,  and 
computation  proceeds  until  there  is  no  contribution.   Hence  timing 
for  the  elementary  functions  will  decrease  somewhat  as  b  is  increased, 
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Checkout 

Several  versions  of  MPP  have  been  used  for  a  total  of  many 
hours  on  IBM  704' s,  7090° s,  and  7094* s. 


Description  of  the  Pseudo-Instructions 

1.  P.I.  will  be  used  to  mean  a  pseudo-instruction  interpreted 
by  this  package. 

2.  M.P.N,  will  be  used  to  mean  "multiple  precision  number." 

3.  The  address  of  the  M„P.N.  in  cells  A,  A  +  1,  A  +  2  shall 
be  A,  and  (A)  specifies  the  M.P.N,  at  address  A. 

4.  A  three  cell  pseudo-accumulator  will  be  designated  by  ACC 
and  (ACC)  refers  to  the  M.P.N,  in  the  ACC. 

5.  ACC1  will  refer  to  the  first  cell  of  the  ACC. 

6.  B  will  indicate  a  single  location  and  (B)  may  be  any  36 
bit  configuration. 

7.  (ACCl)f  means  (ACC1)  treated  as  a  machine  floating  point 
number . 

8.  (ACCl)i  means  (ACCl)  treated  as  an  integer  i  and 
o  <  |i|  <    235. 

9.  (j)  4—  (k)  means  that  the  contents  of  j  are  replaced  by  the 
contents  of  k. 

10.   Each  P.I.  is  written  in  the  following  manner: 

i)   The  address  (bits  21-35)  part  of  the  P.I.  contains  the 
pseudo-operation . 

ii)   The  tag  (bits  18-20)  part  of  the  P.I.  contains  a  tag  of 
0,  1,  or  2.   3,  5,  6,  and  7  may  be  used  on  the  7094. 

iii)   The  decrement  (bits  3-17)  part  of  the  P.I.  contains  the 
address . 


iv)   The  prefix  (sign,  1,  and  2  bits)  is  always 


zero. 
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11.  If  a  P.I.  contains  I  in  the  tag  position,  that  P.I.  is 
indexable  with  either  a  zero,  one,  or  two.   Index  4  must 
not  be  used  while  in  the  interpretive  mode.   3,  5,  6,  and 
7  may  be  used  on  the  7094. 

12.  If  nothing  is  written  in  tag  and  decrement  parts  of  a  P.I., 
then  that  part  of  the  P.I.  may  be  utilized  by  the  programmer. 


List  of  the  Pseudo-Instructions 

In  the  following  list,  the  general  form  of  each  P.I.  is  followed 
by  its  function  and  restrictions.   If  not  specifically  mentioned, 
the  next  P.I.  to  be  interpreted  is  the  next  one  in  sequence. 

:  (ACC)-e  "(A). 

:  (ACC)<_  (A)  . 

:  (A)  <-  (ACC)  . 

:  (ACC)«-  (ACC)+(A). 

:  (ACC)<-  (ACC) -(A) . 

:  (ACC)<—  (ACC) * (A) . 

:  (ACC)4r—   (ACC)/(A),  error  message  if  (A)=o. 

:  (ACC)< (ACC)  . 

:  (ACC)<—  |  (ACC)  |  . 

:  (ACC)« |  (ACC)  |  . 

:  Next  P.I.  is  taken  from  location  B  if 
(ACC)  =  ±  o. 

:  Next  P.I.  is  taken  from  location  B  if 
(ACC)>  +  o. 

:  Next  P.I.  is  taken  from  location  B  if 
(ACC)<   -  o. 

:  Next  P.I.  is  taken  from  location  B  if 
(ACC)  ^  ±  o. 


1. 

CLS , I , A 

2. 

CLA,I,A 

3. 

ST0,I,A 

4. 

ADD ,  I ,  A 

5. 

SUB ,  I ,  A 

6. 

MPY,I,A 

7. 

DVH,I,A 

8. 

CHS 

9. 

SSP 

10. 

SSM 

11. 

TZE , o , B 

12. 

TPL,o,B 

13. 

TMI  ,  o  ,  B 

14. 

TNZ , o , B 

SI- 


15.   ETRA,I,B 


16.   ETZE,I,B 


Departure  is  made  from  the  interpretive 
mode  and  control  is  sent  to  location  B. 

Same  as  15  if  (ACC)=  -  o. 


17.   ETMI,I,B 


18.   CLS1,I,B 


19.   CLA1,I,B 


20.   ST01,I,B 


21.   DVH2N,o,N 


22.   MPY2N,o,N 


23.   STZ,I,A 


24.   AT0B,I,B 


25.   DPEXIT 


26.   DPFLTI 


27 .   DPFLTF 


28.   DPTOM 


29.   DPSQRT 


30.   DP EXP 


:  Same  as  15  if  (ACC)<  -  o„ 

s  (ACC1)<-  -(B) . 

s  (ACC1)<=  (E) . 

:  (B)«=  (ACC1)  . 

s  {ACC)<—  (ACC)/2N   N>o. 

:  (ACC)«-=  (ACC)-2N,  N>o. 

:  (A)'-#=  o„ 

i  (ACC)«™  (ACC)**(B),   j(B)|<  235, 

s  Departure  is  made  from  the  interpretive 
mode  and  control  passes  to  the  next  in- 
struction in  sequence. 

:  (ACC)  is  replaced  by  the  M.P.N,  equivalent 
of  (ACCl)i. 

:  (ACC)  is  replaced  by  the  M.P.N,  equivalent 
of  (ACCl)f. 

t    (ACC1)  is  replaced  by  the  machine  floating 

equivalent  of  (ACC) .  Numbers  out  of  machine 

range  are  replaced  by  extreme  values,  i.e., 
o  or  (377777777777)8- 

:  (ACC)  is  replaced  by  the  square  root  of 
|  (ACC)  j  . 

:  (ACC)  is  replaced  by  the  exponential  of  (ACC) 
The  argument  must  be  less  than  2°~*-    ln2  in 
absolute  value.   For  b=9  the  absolute  value 
of  the  argument  must  be  less  than  177.4. 
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DPLNX 


32.  DPCOS 

33.  DPSIN 

34.  DPATN 

35.  DPASIN 

36 .  DPACOS 

37.  DPSINH 

38.  DP COSH 

39.  LDQ,I,A 

40.  STQ,I,A 

41 .  CACMQ 

42.  DPFBCD,I,A 


(ACC)  is  replaced  by  the  natural  logarithm 
of  (ACC) .  If  the  argument  is  less  than  or 
equal  to  zero,  no  computation  is  done,  and 
the  (ACC)  is  unchanged. 

(ACC)  is  replaced  by  the  cosine  of  (ACC) . 
The  argument  plus  ff/2  must  be  less  than  2^5 
in  absolute  value. 

(ACC)  is  replaced  by  the  sine  of  (ACC) . 
The  absolute  value  of  the  argument  must 
be  less  than  235 

(ACC)  is  replaced  by  the  arc  tangent  of  (ACC) . 
The  result  is  between  -  TT^  and  ir/2. 

(ACC)  is  replaced  by  the  arc  sine  of  (ACC) . 
The  absolute  value  of  the  argument  must  not 
be  greater  than  1. 

(ACC)  is  replaced  by  the  arc  cosine  of 
(ACC) .   The  absolute  value  of  the  argument 
must  not  be  greater  than  1. 

(ACC)<—  hyperbolic  sine  of  (ACC) . 

(ACC)  <«—  hyperbolic  cosine  of  (ACC). 

(Pseudo-MQ) <—  (A)  . 


(A)< 


(pseudo-MQ) . 


43 .  GET3 , I , R 

44 .  PUT3 , I , R 


Next  P.I.  taken  if  (ACC)>  (pseudo-MQ). 
Skip  one  P.I.  if  (ACC)  =  (pseudo-MQ).   Skip 
two  P.I.s  if  (ACC) <  (pseudo-MQ) . 

(ACC)  is  replaced  by  the  M.P.N,  equivalent 
of  the  BCD  number  stored  starting  at  loca- 
tion A.   This  P.I.  provides  a  way  of  entering 
constants  into  a  code.   An  example  is  given 
in  the  section  on  input. 

(ACC)«-  (R) . 

(R)<-  (ACC) . 
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45.  GET1,I,R       :  (ACC1)«—  (R)  . 

46.  PUT1,I,R       :  (R)  <—  (ACCl)  . 

47.  TXI.D,B       :  Index  register  1  is  incremented  by  P  and 

the  next  P.I.  is  taken  from  location  B. 

Pseudo-instructions  43,  44,  45,  and  46  above  provide  a  crude 
way  of  communicating  with  MPP  when  the  package  is  assembled  as  a 
relocatable  subprogram.   R  is  the  address,  relative  to  zero,  of 
some  1  or  3  cell  value  within  MPP. 

Pseudo-instructions  48  through  55  listed  below  have  proven  of 
little  value  but  are  included  for  the  benefit  of  codes  already 
written.   The  function  of  each  P.I.  is  described  by  giving  the 
machine  instruction  which  the  P.I.  simulates.   K=0,  1,  or  2.   K 
may  also  be  3,  5,  6,  or  7  on  the  7094. 

48.  TIX,D,B  :  TIX  B,1,D 

49.  TXH,D,B  :  TXH  B,1,D 

50.  TNX,D,B  :  TNX  B,2,D 

51.  TXL,D,B  :  TXL  B,2,D 

52.  LXD,K,B  :  LXD  B,K 

53.  SXD,K,B  :  SXD  B,K 

54.  LXA,K,B  :  LXA  B,K 

55.  SXA,K,B  :  SXA  B,K 

Example ;   Given  the  argument  x  in  the  ACC,  compute 


-16xe  X/4   Varc  tan(lnx2) 
and  leave  the  result  in  the  ACC.   Assume  )x|  >  1. 
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YFUN  SXD  YFUN+4,4 
TSX  ENTRY ,4 
STO,0,X 
MPY,0,X 
DPLNX 
DPATN 
DPSQRT 
STO,0,T 
CLS,0,X 
DVH2N,0,2 
DPEXP 
MPY,0,T 
MPY,0,X 
MPY2N,0,4 
CHS 

DP EXIT 
LXD  YFUN+4,4 
FINIS  TXI  ENTRY, 4,-1 
X  BSS  3 
T  BSS  3 


SAVE  INDEX  4 
ENTER  INTERPRETIVE 
SAVE  ARGUMENT 
SQUARE  ARGUMENT 
LN  OF  X**2 
ARC  TANGENT 
SQUARE  ROOT 
SAVE  PART  OF  RESULT 

-X/4 
EXPONENTIAL 


LEAVE  INTERPRETIVE 


GO  GET  NEXT  P.I. 


The  above  coding  could  be  executed  with  the  P.I 


YFUN, 


Method  and  Accuracy  For  Elementary  Function  Pseudo-Instructions 
The  argument  (ACC)  -  x  =  F-2^  =  g-lO^  ,  .1<  |  g  |  <  1 
1.   Square  Root; 

An  initial  approximation  yQ  is  computed  as  follows: 
Case  1:   p  even,  y0  =  %(  |  f|  +1  )-2p/2 

£+1 

Case  2:   p  odd,  yQ  =  H(     |  F  |  +  h    )-2   7~ 
The  square  root  is  then  yc  obtained  from  the  recursion  formula 

n+1      n    yn 
The  result  should  be  correct  to  .3(105-b)  significant  decimal  digits 
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2.  Exponential; 

Put  x/ln2  =  I  +  f,  I  an  integer  and  o  <  |  f |  <  1. 

_     x    „I    f  ln2 
Then  e   =  2  •  e 

e      is  computed  from  the  exponential  series.   The  result 
should  be  good  to  .3(105-b)-l  significant  decimal  digits  if  |x|  <    1. 
If  |x|  >  1,  the  result  should  be  correct  to  .3 (105-b) -(q+1)  signifi- 
cant decimal  digits. 

3.  Natural  Logarithm 
In  x  =  In  F  +  p  ln2. 

in  P  =  2  (<f=i)  ♦  l/3(f^)3  +  l/5(f^)5  ♦  ...}  . 

The  result  should  be  correct  to   .3 (105-b) -1  significant  decimal 
digits,  except  when  a  loss  of  significance  occurs  in  computing  F-l. 

4.  Sine  and  Cosine 

Put  x/if/4   =  I    +   f,    I    an  integer   and   o  <    |   f  |    <    1. 

sin  x   =  sin    (I'T  +   f'^/4). 

Either   sin   f"ff/4  or  cos    f*if/4   is   computed    from  the   sine 
or   cosine   series . 

Cos  x  is  obtained  by  computing  sin  (x  +  -If/2)  . 

The  accuracy  of  sin  x  and  cos  x  should  be  the  same  as  that 
given  for  the  exponential  function,  except  for  cos  x  when  a  loss 
of  significance  occurs  in  computing  x  +  TT/2. 

5.  Arc  Tangent 

The  following  relations  are  used, 
arc  tan  x 


^a-         2n+l 
=  y        (-Dn  X2n+1  *■■—  -  arc  tan  (-)=  -  arc  tan  (-x) 


n=o 

x-a 


=  arc  tan  a  +  arc  tan 


1+ax  ' 


The  result  should  be  correct  to  .3 (105-b) -1  significant 
decimal  digits. 
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6.  Arc  Cosine 

arc  cosine  x  =  arc  tangent  (  y  1-x  /x)  ,  l>|x|>o;  arc 
cosine  o  ■  TT/2 . 

The  result  should  be  correct  to   .3(105-b)-l  significant 
decimal  digits,  except  when  significance  is  initially  lost  in  com- 
puting 1-x  . 

7 .  Arc  Sine 

arc  sine  x  =  1T/2  -arc  cosine  x. 

8.  Hyperbolic  Cosine  and  Hyperbolic  Sine 
cosh  x  =  (e   +  e   )  /2. 

sinh  x  =  (e   -  e   )  /2. 

Results  should  have  the  same  accuracy  as  that  described 
for  the  exponential  function,  except  for  sinh  x  when  x  is  close  to 
zero  and  loss  occurs  in  subtraction. 

9 .  x*9     with  the  P . I .  :   ATOB ,  I ,  B 

If  (B)=/0has  n  significant  binary  bits  and  k  of  these  are 
non-zero  than  x*  is  computed  using  n+k-2  multiplies.  One  division 
is  required  if  /3  is  negative.  This  P.I.  is  used  only  for  raising 
a  multiple  precision  number  to  a  integral  power.   o^so,  X°  =  l« 

Note;  It  is  the  user's  responsibility  to  insure  that  the 
argument  x  is  in  the  proper  range  before  computing  Vx~,  lnx,  sin  x, 
cos  x,  sin  x,  cos~1x. 


Significance  and  Number  Range  Table 

Let  a  multiple  precision  number  have  the  following  representa- 
tion: 

N=f-2P  =  F-10q  ,    h   <    I  f  I  <  1  ,  or  f  =  o, 

.1   <     |F|  <  1  ,  or  F  =  o. 
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The  following  table  gives  a  few  values  of  b,  the  approximate 


upper  limits  for  |p  and  q  and  number  of  significant  decimal 
digits  in  the  number.  As  described  previously,  b  is  the  number  of 
bits  allocated  for  the  binary  power  of  two  of  a  M.P.N. 


Number  of 

1  P 

must  be 

1 q   must  be 

Significant 

b 

less  than 

less  than 

Decimal  Digits 

9 

256 

78 

28.8 

10 

512 

155 

28,5 

11 

1024 

309 

28.2 

12 

2048, 

617 

27.9 

13 

4096 

1234 

27.6 

14 

8192 

2467 

27.3 

15 

16384 

4933 

27.0 

16 

32768 

9865 

26.7 

17 

65536 

19729 

26.4 

A  few  octal  representations  of  numbers  when  b 

000000000000 
400000000000 
314631463146 


=  9  are: 


1  =  h    '    2  =     200000000000 

■5=-(ij+l/8) -23=  640000000000 
-1 


3=(.6) -2 


231463146314 


000000000401 


400000000403 


146314632377 


Errors  Detected 

There  are  6  types  of  serious  errors  detected  in  the  package. 

El:  An  error  condition  exists  if  the  unnormalized  fractional 
part  of  a  product  is  less  than h  in  absolute  value.  The 
error  cannot  occur  if  only  normalized  M.P.N.s  are  used. 
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E2 :   Overflow  which  may  occur  in  several  parts  of  the  package 
or  in  the  program  of  the  user.   Use  of  the  MODIFY  sub- 
routine will  usually  correct  overflow  problems.   Results 
which  underflow  are  automatically  set  to  zero. 

E3 :   Division  by  zero  is  not  allowed. 

E4:   Division  by  an  unnormalized  M.P.N,  is  not  allowed. 

E5 :   An  error  condition  exists  if  an  end  of  file  is  read  by 

the  input  subroutine.   The  input  tape  will  be  repositioned 
in  front  of  the  end  of  file  before  taking  the  error  exit. 

E6:   Some  format  errors  are  detected  during  output. 

If  one  of  the  above  errors  is  detected,  a  message  is  written 
on  the  output  tape  and  the  instruction  STR  is  executed.   The  STR 
will  normally  cause  control  to  go  to  the  dump  routine  of  IBSYS. 

Error  Exit  Modification 

Any  or  all  of  the  above  error  exits  may  be  modified  by  use 
of  the  sequence: 


<X, 

TSX 

SETERS ,  4 

<*+l 

PZE 

EXIT1 

X+2 

PZE 

EXIT2 

*+3 

PZE 

EXIT3 

<*+4 

PZE 

EXIT4 

*+5 

PZE 

EXIT5 

«+6 

PZE 

EXIT6 

*+7 

PZE 

NOBSR 

<x+8 

TRA 

CEXIT 

<*+9 

Return 

This  sequence  sets  error  Ei  to  transfer  control  to  EXITi  and 
the  error  message  is  not  printed.   Ei  is  unchanged  if  EXITi=*0.   If 
NOBSR^O  then  the  tape  is  not  backed  over  the  end  of  file  during  E5 , 
If  the  contents  of  °<+8  are  not  zero  then  the  STR  exit  is  replaced 
by  the  contents  of<*+8. 

Under  the  FORTRAN  Monitor  under  IBSYS,  it  may  be  preferable 
to  have  (<*+8)=  CALL  EXIT  in  order  to  avoid  returning  to  IBSYS  and 
dumping . 
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Subroutines  in  the  Package 

1.  Number  Range  Modification  Subroutine  -  MODIFY 

If  it  is  necessary  to  have  a  greater  range  of  numbers  than 
is  possible  when  b  =  9,  then  b  can  be  increased  by  giving  the  calling 
sequence 

*     TSX  MODIFY, 4 
*(  +1   PZE  B 
oi  +2   return 

This  subroutine  sets  up  the  entire  package  so  that  future  com- 
putations will  have  (105-B)  bits  for  the  fractional  part  of  numbers 
and  B  bits  are  used  for  the  power  of  2  exponent. 

There  may  be  constants  outside  this  package  which  have  to  be 
modified  also.   This  can  be  accomplished  by  inserting  a  transfer 
into  MODIFY  +  48  before  giving  the  above  calling  sequence.   This 
transfer  should  be  to  a  code  which  will,  for  each  M.P.N,  to  be  mod- 
ified, clear  and  add  logical  the  third  word  of  the  M.P.N,  and 
TSX  ADJUST, 4.   Control  will  return  to  the  location  following  the 
TSX,  at  which  point  the  accumulator  should  be  stored  logically  back 
into  the  third  cell  of  the  M.P.N. 

2.  Subroutine  for  Solution  of  Linear  Equations  -  EQS 

This  subroutine  solves  a  system  of  N  simultaneous  linear 
equations  by  the  method  of  elimination  with  pivot  selection. 

To  solve  AS=B  for  S  the  following  calling  sequence  is  given. 

o(  TSX  EQS,  4 

°L+1  PZE  A,0,S 

o(+2  PZE  N  (N=  number  of  equations) 

°C  +3      error  return  (singular  matrix) 

o(  +4      normal  return 

A,  B,  and  S  are  stored  as  follows s 

aj_j  is  stored  in  location  A+3  ( j-l+(i-l)  (N+l)). 

bk  is  stored  in  location  A+3 (k (N+l) -1) ,  k=l,2,...,N. 

sm  will  be  stored  in  location  S+3(m-l),  m=l,2,...,N. 

The  matrix  A  is  destroyed  during  the  computation. 

It  should  take  about  2  minutes  to  solve  50  equations  on  the  7090. 
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Input  Subroutine  -  DPLOAD 


This  subroutine  reads  and  converts  BCD  information  from 
SYSINl.   Records  consisting  of  DEC, OCT, BCD, BCI ,  or  TRA  type  infor- 
mation are  allowed.   In  addition,  if  column  1  is  not  punched  with 
D  and  columns  8,  9,  and  10  are  blank  then  DEC  is  understood. 

M.P.N,  records  are  also  handled  by  this  subroutine,  and  they 
are  identified  by  having  a  D  in  column  1  followed  by  the  multiple 
precision  numbers  starting  in  column  2.   Multiple  precision  numbers 
are  written  exactly  as  numbers  on  DEC  records,  except  up  to  30 
significant  decimal  digits  may  be  written  for  each  number;  and  the 
absolute  value  of  the  decimal  power  of  ten  exponent  (n)  can  be  any 
integer  less  than  one  third  the  core  size  minus  the  number  of  digits 
to  the  right  of  the  decimal  point  in  the  fractional  part  of  the 
number.   A  table  of  conversion  constants  for  I  n  |  ^  68  is  provided. 
If  |  n J  >  68  the  conversion  constant  is  computed,  requiring  extra 
conversion  time  and  may  introduce  round  off  errors. 

All  integers  are  converted  exact  by  this  subroutine,  while  non- 
integers  are  correct  (unrounded)  to  the  full  significance  in  effect 
at  time  of  conversion. 

Loading  is  accomplished  by  giving  the  calling  sequence 

C>£      TSX  DPLOAD,  4 

o^+l    P  fi.O.V 

o(  +2        return 

where  P  is  either  PZE  or  TXL. 

OCT,  BCD,  DEC,  BCI,  or  D  type  records  are  read  and  converted 
consecutively  until  the  storage  is  filled  from  /@   through  )f   .   It 
must  be  remembered  that  multiple  precision  numbers  take  3  cells  while 
all  other  types  take  one  cell  each.   A  TRA  ^,1  record  will  terminate 
loading  and  control  will  be  transfered  as  specified  by  the  TRA.   Also 
a  TOC  *L,I  is  recognized  while  loading  and  has  the  same  function  as 
TRA  except  that  L  is  treated  as  an  octal  integer.   Records  may  be 
up  to  120  characters  in  length. 

DECOCT,  BCI,  and  BCD  records  may  have  a  decimal  address,  or 
the  character  equals  (=)  followed  by  an  octal  address,  in  columns 
1  through  6  of  the  record.   In  either" case  the  loading  of  data  begins 
at  the  new  address. 
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DPLOAD  checks  for  a  few  clerical  errors  which  may  be  in  the 
data  field  of  DEC  and  D  type  records.   If  P  of  the  calling  sequence 
is  zero  the  machine  stops  after  detecting  an  error.   If  P  =  7,  the 
machine  does  not  stop  on  error  detection. 

Each  error  detected  in  a  record  causes  the  entire  record  to  be 
printed  on  line  with  the  column  number  in  error  and  the  number  of 
records  which  have  been  read  up  to  this  point.   If  the  machine  stops 
(i.e.,  P  =  o)  after  an  error  print  it  may  be  restarted  and  loading 
will  continue  as  if  P  =  7. 

After  reading  a  block  of  data  with  DPLOAD,  a  check  may  be  made 
to  determine  if  any  errors  were  detected.   This  is  accomplished  by 
giving  the  instruction 

°L  TSX  TEST,  4 

Control  returns  to   o(  +1  with  the  accumulator  address  containing 
the  number  of  records  read  by  DPLOAD  since  the  last  time  TEST  was 
entered.   The  accumulator  sign  is  positive  if  no  errors  were  de- 
tected, negative  if  one  or  more  errors  were  detected. 

The  errors  detected,  and  actions  taken  if  the  program  is  con- 
tinued are  listed  below. 

Error  Action  Taken 

1.  Illegal  character  (i.e., 
anything  other  than 

o ,  1 ,..,.,  9  ,  E ,.,+,- ,  dash ,  B , 

or  comma)  ignored 

2.  2  consecutive  decimal 

points  2nd  point  ignored 

3.  2  consecutive  commas       zero  separates  the  commas 

4.  point  followed  by  plus     plus  followed  by  point 

5.  point  followed  by  minus    minus  followed  by  point 

6.  point  followed  by  dash 

(dash  is  treated  as  minus)  dash  followed  by  point 

7.  comma  followed  by  blank    comma  followed  by  zero, 

followed  by  blank 
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(continued) 
Error 


Action  Taken 


8. 


an  operation  in  columns 
8,  9,  10  other  than 
DEC;OCT/BCD/TRA,BCI,  or 
TOC 


entire  record  ignored 


Example 


Suppose  b=9,  and  the  following  calling  sequence  was  given  starting 

at  (50060)_  =  (20528)... 

8  10 

TSX  DP LOAD, 4 
TXL  INPUT, 0,-1 

The  14  cards  listed  below  are  to  be  read  from  SYSINl. 

BCI    1, +40000 
BCD  1000000 

257X 
DEC  .  .-1, ,240 
BDC  II DENT 
D-.314159265358979323846264339C28E01 
DO. 

DEC  2. ,0, 
448  OCT  460000000355 
=611  OCT  450100000437,060100000623 
TOC  50060 
TRA  20528 
TOC  0,4 
TRA  2,4 

An  on-line  print  out,  except  for  extra  spacing,  would  appear  as 
follows . 


DEC  257X 
DEC  ..-1, ,240 
DEC  ..-1,,240 
DEC  ..-1, ,240 
BDC  II DENT 
D-.314159265358979323846264339C28E01 
DEC  2 . , 0 , 


CHECK  CARD  3  IN  COLUMN  15 
CHECK  CARD  4  IN  COLUMN  13 
CHECK  CARD  4  IN  COLUMN  14 


CHECK  CARD  4 
CHECK  CARD  5 
CHECK  CARD  5 
CHECK  CARD  7 


IN  COLUMN  17 
IN  COLUMN  8 
IN  COLUMN  31 
IN  COLUMN  17 
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core  storage,  starting  at  symbolic  location  INPUT,  would  contain  the 

following  octal  words. 

200400000000  000000000000  000000000401  575631463146  000000000000 
000000000360  711037552421  413214106461  461461510402  000000000000 
000000000000  000000000000  202400000000  000000000000  000000000000 

The  3  octal  words  on  the  ninth  and  tenth  cards  would  be  stored  in 

locations  (488).,,,  (611) _  and  (612) _,  and  control  would  go  to  octal 

10        8  8 

location  50062. 

Note  that  when  an  entire  record  is  rejected  the  record  (card)  count 
is  not  incremented. 

The  variable  field  for  BCI,OCT,DEC,TRA,  and  TOC  type  records 
may  begin  anywhere  in  columns  12  though  16. 

Input  by  DPLOAD  is  accomplished  by  use  of  the  subroutine  RTX 
(described  below)  which  in  turn  uses  IOEX.   All  input  is  single 
record  buffered  so  it  is  recommended  that  the  user ' s  data  deck  have 
a  dummy  data  card  at  the  end.   However,  an  end  of  file  will  usually 

be  sufficient. 

DPFBCD 

DPLOAD  is  used  for  execution  of  the  pseudo  instruction 
DPFBCD, I,A.   The  conversion  is  accomplished  by  forcing  DPLOAD  to 
get  its  input  from  storage  rather  than  SYSIN1. 

Only  one  number  can  be  converted  by  DPFBCD  by  a  single  use. 
The  BCD  number  is  written  as  D  type  data  and  the  number  mugt  be 
preceded  by  the  character  D  and  terminated  by  a  blank. 

Constants  can  be  entered  into  the  code  in  BCD  form  and  con- 
verted in  this  manner  as  an  initial  operation  following  any  use  of 
the  MODIFY  subroutine, 

Example  for  DPFBCD 

The  following  code  could  be  used  to  enter  the  constants 
1.,-tt/4,  and  ln2  into  storage  locations  ONE,  MPIOV4,  and  L0G2. 
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TSX     ENTRY,  4 

DPFBCD,0,C2 

STO,0,LOG2 

DPFBCD,0,C1 

STO,0,ONE 

DPFBCD,0,C3 

STO,0,MP1CP74 

DP EXIT 


CI      BCI 

2,D 

C2      BCI 

6,D 

C3      BCI 

6,D 

LOG 2    BSS 

3 

MPIOV4  BSS 

3 

ONE     BSS 

3 

Readinq 

Errors 

D10E-1 

06931471805599453094172321E1 
D-. 7853981633974483096156608458 


Noise  records  (records  of  length  less  than  3  words)  are  ignored 
IOEX  prints  a  message  if  a  record  is  permanently  redundant  and  the 
record  as  read  the  last  time  is  accepted.   Reading  through  an  end 
of  file  is  considered  an  error,  but  this  can  be  changed  by  use  of 
the  SETERS  routine. 


4.   Input  Record  Scan  Deletion  Subroutine  -  NOSCAN 

The  DPLOAD  routine  can  be  set  to  omit  error  prints  1. 
through  8.  by  giving  a  TSX  NOSCAN, 4  but  the  actions  taken  on  these 
errors  are  the  same  as  described  before. 
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5.   BCD  Read  Subroutine  -  RTX 

This  subroutine  reads  BCD  input  records  from  SYSINl. 
The  calling  sequence  for  RTX  is  as  follows. 

of.  TSX     RTX,  4 

ot+1      PZE     A,0,EOFR 
cC  +2      Return 

One  BCD  record  or  the  first  22  words  of  a  BCD  record  is  read 
into  core  storage  starting  at  symbolic  location  A. 

Records  are  single  record  buffered  automatically.   An  end  of 
file  causes  control  to  be  sent  to  symbolic  location  EOFR.   If  EOFR 
were  *-l  or  0  then  end  of  files  would  be  ignored  with  buffering 
initialized  at  the  beginning  of  each  file  to  be  read. 

RTX  uses  IOEX  for  reading.   Noise  records  are  ignored,   Per- 
manent redundancy  causes  IBSYS  to  print  a  message  and  then  RTX 
will  accept  the  record  as  read  the  last  time. 


Warning 

Special  care  is  required  when  non-IOEX  input/output  routines 
are  used  in  conjunction  with  MPP  input/output.   The  FORTRAN  II 
library  routines  are  examples  of  a  type  which  will  cause  conflict, 
Such  routines  may  be  used  only  after  the  channels  have  been  freed 
by  IOEX.   The  sequence 

TCOX  * 
TCOX  *-l 

for  each  channel  X  is  sufficient. 
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6.   Binary  to  Decimal  Conversion  and  Output 

a.   Output  Subroutine  -  DPBDC3 

This  subroutine  converts  information  in  core  and  prints 
it  on  SYSOU1.   Binary  integers  may  be  printed  in  octal  or  decimal 
form.   Multiple  precision  numbers  may  be  printed  in  floating  or  fixed 
point  decimal  form.   Integers  are  exact  while  multiple  precision  num- 
bers are  rounded  to  a  specified  number  of  places.   BCD  information 
can  be  printed  also. 

Output  is  accomplished  by  giving  the  calling  sequence 

oi  TSX    L,4 

a(+l         PZE     A,0,B 
oi  +2   Return 


L      TRA     DPBDC3 

BCI     n , format 

The  block  of  numbers  in  continuous  core  storage  from  A  through 
B  will  be  printed  as  specified  by  the  format  and  control  returns  to 
ci   +2.   L  is  any  core  location  and  the  BCD  format  (see  format  descrip- 
tion) may  be  of  any  length. 


Example :   The  following  code  would  cause  100  full  word  integers 
stored  starting  at  location  Y  to  be  printed  as  decimal  integers 

Q    TSX   F,4 

Q+l   PZE    Y,o,Y+99 

Q+2  Return 


TRA    DPBDC3 
BCI    1,8114 
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b.   Output  Subroutine  -  D(BDC) 

This  subroutine  handles  conversion  in  exactly  the  same 
way  as  DPBDC3 .   The  only  difference  is  in  the  way  information  is 
picked  up  from  core  for  conversion  and  printing.   D(BDC)  allows 
printing  of  information  which  is  not  stored  sequentially  in  core. 

Output  is  accomplished  by  giving  the  calling  sequence 

o^      TSX    L,4 

o(  +1    PRF     OP,I,A 


<^  +k    PRF     OP  ,  I ,  A 
TRA     D(FIL) 


L      TRA     D(BDC) 

BCI    n , format 

The  user's  coding  in  <<+l  through  o(+k  determines  which  numbers 
are  to  be  converted  and  printed  as  specified  by  the  format.   See 
format  description. 

The  block  of  coding  in  o(+l   through  <tf+k  may  not  contain  any 
instructions  other  than  the  following  types. 

i.   Any  machine  instruction  which  does  not  alter 

sequential  flow  of  control.   For  example,  test  and 
unconditional  transfer  instructions  are  not  allowed. 

ii.   Any  of  the  instructions  TXI,  TIX,  TXH,  TNX,  TXL. 

iii.   The  STR  instruction,  which  has  the  pseudo  operation 
CLA,I,A  in  its  address,  tag,  and  decrement  portions. 
This  is  the  instruction  used  to  pick  up  numbers  for 
conversion  and  printing. 

Index  register  4  may  not  be  used  in  the  above  instructions. 
All  other  index  registers  are  available,  but  only  1  and  2  are 
saved  and  restored. 
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Example 


The  following  code  would  cause  the  single  celled  number  N 
and  the  3  celled  numbers  Bi ,  i  =  1, 2 , . . . , (N)  to  be  printed  according 
to  some  format. 


TSX 

FORM, 4 

STR 

CLA,0,N 

hXA 

N,l 

AXT 

0,2 

STR 

CLA,2,B 

TXI 

*+l,2,-3 

TIX 

*-2,l,l 

TRA 

D(FIL) 

FORM 


TRA 
BCI 


D(BDC) 
n , format 


N 
B 


BSS 
BSS 


1 
30 


Note 

The  store  and  trap  (STR)  instruction  used  in  the  calling 
sequence  is  not  executed,  but  is  used  to  isolate  types  i  and  ii 
instructions  from  the  data  pickup  instruction.   All  instructions 
between  the  TSX  and  TRA  D(FIL)  are  handled  interpretively. 


Warning 

Special  care  is  required  when  non-IOEX  input/output  routines 
are  used  in  conjunction  with  MPP  input/output.   The  FORTRAN  II 
library  routines  are  examples  of  a  type  which  will  cause  conflict, 
Such  routines  may  be  used  only  after  the  channels  have  been  freed 
by  IOEX.   The  sequence 

TCOX  * 
TCOX  *-l 

for  each  channel  X  is  sufficient. 
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Format  Description 

A  format  consists  of  any  number  of  field  specifications. 
The  general  field  specification  is  of  the  following  form. 

o(S  /3  (tf)gPnTd-wZ 

ei.i/3  i   Jf  ,g/n,d,  and  w  are  positive  integers. 

T  is  the  type  conversion  character  and  T  is  I,  E,  F,  0,A, 
or  X. 

Z  is  a  slash,  comma,  or  blank,  and  Z serves  as  a   field 
specification  separator. 

Each  non-numeric  Hollerith  character  in  a  field  specification 
indicates  what  is  to  be  done  with  the  integer  just  preceding  it. 

Each  field  specification  is  executed  n  times  and  n  data  fields 
are  printed,  each  having  a  total  column  width  w  with  the  spread  of 
information  being  determined  by  o(,  fl  ,Y,g,  and  d. 

T  =  I  specifies  full  word  binary  to  decimal  integer  conversion. 

T  =  E  specifies  binary  to  floating  point  decimal  conversion 
of  multiple  precision  numbers. 

T  =  F  specifies  binary  to  fixed  point  decimal  conversion  of 
multiple  precision  numbers . 

T  =  0  specifies  full  word  binary  to  octal  conversion. 

T  =  A  specifies  BCD  to  BCD  conversion. 

T  =  X  specifies  insertion  of  blank  characters. 

If  either  of  the  integers  o(, /8  ,?{,    or  g  is  zero  then  that 
integer  and  the  following  Hollerith  character  may  be  omitted. 

The  format  is  scanned  from  left  to  right.   Conversion  as 
indicated  by  a  specification  is  completed  before  checking  on  the 
separator  Z.   If  Z  is  a  blank,  then  all  information  converted  at 
this  point  is  printed  and  if  there  is  still  more  information  to 
be  converted  and  printed,  scanning  goes  back  to  the  beginning  of 
the  format.   If  Z  is  a  slash,  all  information  converted  at  this 
point  is  printed  and  scanning  for  the  next  line  starts  immediately 
after  the  slash.   K  consecutive  slashes  separating  two  specifications 
cause  K-l  blank  lines  to  be  printed.   If  Z  is  a  comma,  conversion 
continues  with  the  specification  following  the  comma  being  scanned. 
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Binary  to  Decimal  Integer  Conversion  (T=I)   i  gPnlwZ  7 

n  single  celled  binary  integers  are  converted  and  printed  as 
decimal  integers,  each  with  a  total  column  width  w;  and  each  number 
is  enclosed  in  parenthesis  if  g^o.   The  restrictions  are  n>o,  wi  26 

Example :   Suppose  the  binary  equivalent  of  the  integers  1,  -2,  3, 
7000,  and  -56789  are  in  consecutive  storage  starting  at  location  C. 
Then  the  sequence: 

TSX  FORM, 4 

C,o,C+4 


FORM  TRA  DPBDC3 

BCI  2,1P1I4,4I5 

would  cause  the  following  line  to  be  printed. 

(1)    -2     3  700056789 

Note  that  the  negative  sign  of  the  last  number  would  not  be  printed 
since  w  —  5. 

Binary  to  Floating  Point  Decimal  Conversion  (T=E)  ^<<S  Q  (y)gPnEd«wZ  I 

n  3  cell  multiple  precision  numbers  are  converted  and  printed, 
each  number  having  a  total  column  width  w.   Information  within  the 
column  is  arranged  as  follows  from  left  to  right. 

1.  w-«<-tf-V-d-2    blanks. 

2.  q,  the  decimal  exponent  of  the  number  occupying  /Q   places 
and  in  parenthesis  if  g^o. 

3.  The  sign  of  the  number  (blank  for  +,  or  a  minus  sign). 

4.  d.     significant  digits  of  the  number. 

5.  A  decimal  point. 
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Res 


6.  d  significant  decimal  digits  (rounded). 

7.  q,  the  decimal  exponent  of  the  number  occupying  y  places 
and  in  parenthesis  if  g^o. 

frictions  are  n  >■  o,  w>o^+/<£+y+d  +  2,  fi*  26,  y^  26,  d  +  oi<30 


Example:   The  format  1S4(1P1E5 . 13 , 3) 1E10. 19  would  cause  the  binary 
equivalent  of  12345.67  and  -.00765432  to  be  printed  as  follows. 

(4)  1.23457     -.7654320000  -2 

Binary  to  Fixed  Point  Decimal  Conversion  (T=F)       i   «(SnFd.wZ  ? 

n  3  cell  multiple  precision  numbers  are  converted  and  printed 
as  fixed  point  decimal  numbers ,  each  number  having  a  total  column 
width  w.   Information  within  the  column  is  arranged  from  left  to 
right  as  follows.   Before  conversion,  the  number  is  multiplied  by 
10* 

1.  The  integral  part  of  the  number  (with  sign)  occupying  w-d-1 
spaces . 

2.  A  decimal  point. 

3.  d  decimal  places  (rounded)  of  the  fractional  part  of  the 
number . 

Restrictions  are  n>o,  w>d  +  2+  number  of  digits  in  integral  part. 

This  subroutine  arbitrarily  refuses  to  print  fixed  point  multiple 
precision  numbers  whose  absolute  value  is  greater  than  2"*5-l.   if 
such  a  number  is  found,  the  specification  is  treated  as  E  type  rather 
than  F  type  and  the  number  is  printed  in  floating  point  form,  but 
the  power  of  ten  exponent  cannot  be  printed  unless  either  ^(or  y) 
is  written  in  the  specification.    /^and  }f   are  ignored  in  F  type 
conversion,  so  no  harm  is  done  by  specifying  one  of  the  two  just  in 
case  large  numbers  are  printed. 


Example :   The  format   3)3F10.20  would  cause  the  binary  equivalent 
of  -101.01,  1040,  and  .00004778966  to  be  printed  as  follows. 

-101.0100000000       .1000000000  41  .0000477897 
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Binary  to  Octal  Conversion  (T=@,    $  o(  S  /3   ^)gPnOd-wZ  £ 

n*  single  celled  binary  integers  are  converted  and  printed 
as  octal  integers,  with  locations  if  desired.   Each  number  has  a 
total  column  width  w  and  information  within  the  column  is  arranged 
from  left  to  right  as  follows. 

L   w- -°{  ~/3   -  y  -  d  blanks. 

2„   The  location  as  a  decimal  integer  occupying  /3   places  and 
in  parenthesis  if  g  =  1,  3,  9S  or  11 , 

3.  The  location  as  an  octal  integer  occupying  y  places  and 
in  parenthesis  if  g  =  2,  3,  10,  or  11. 

4.  4  blanks. 

5.  The  last  d  octal  digits  of  the  number. 

Restrictions  are  d  <  12 ,  /3   <    26, tf  <  26,  w»<?{  +  /?+S'+d,  n>o. 
If  g  >  8,  then  as  far  as  the  calling  sequence  is  concerned,  the 
number  described  above  is  not  counted  as  having  been  printed;  thus 
allowing  a  single  cell  to  be  printed  both  in  octal  and  decimal  with 
or  without  decimal  and/or  octal  locations?  or  if  d  =  0  allowing  a 
location  to  be  printed  with  either  a  single  celled  number  or  multiple 
precision  number. 

Example s   The  format  1S5 (9P10o. 7 , 116, 2S3) 1012.21, 103.7  would  cause 
the  octal  numbers  -17,314631463146  and  7777777  in  decimal  locations 
100,  101,  and  102  to  be  printed  as  follows, 

(100)    -15     145   314631463146     777 

*   If  n  >  1  then  g  must  be  less  than  8. 


BCD  Conversion  (T*A)   j   nAwZ  I 


n  BCD  words  are  printed  as  n  BCD  fields  with  each  field  having 
a  total  column  width  w.   If  w  >  6,  then  each  field  consists  of  w-6 
blanks  followed  by  the  BCD  word  consisting  of  6  characters.   If 
w  <  6  then  the  field  consists  of  the  left  most  w  characters  of  the 
BCD  word. 

Blank  Insertion  (T=X)   £nXZ^ 

n  blanks  are  inserted  into  the  print  line. 
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BCD  Characters  Within  the  Format 

All  of  the  previously  described  field  specifications  may  be 
preceded  by  iC  followed  by  i  Hollerith  characters.   The  i  characters 
will  be  printed  just  before  the  first  number  printed  by  the  specifi- 
cation which  follows. 

Example :   Given  the  binary  equivalent  of  the  integer  2,  and  multiple 
precision  numbers  .006,  -10   ,  -10   ,  and  106  starting  in  decimal 
core  location  100;  the  calling  sequence 

TSX  FPR,4 

100,0,112 


FPR  TRA  DPBDC3 

BCD  623C   TEMPORIES  STARTING  AT5(8Pl0o.6/ 
BCD  51I4,3)2F4.1o,1P1S5(2E5.16 
would  cause  the  following  two  lines  to  be  printed. 
TEMPORIES  STARTING  AT   100   • 

2      .0060  -.1000  16     (-8)-l. 00000      (6)  1.00000 

BCD  Tape  Write  Subroutine  -  WOT 

This  subroutine  is  used  by  MPP  for  all  of  its  printing  of 
BCD  records  on  SYSOUl .   Output  is  single  record  buffered  and 
IOEX  is  trusted  to  write  the  record  correctly.   If  the  end  of  tape 
is  sensed  a  message  is  printed  and  the  machine  pauses  for  a  fresh 
SYSOUl. 
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Printing  of  Tables        -fmKc<  S/3(y)gPnTd-wZj- 

To  make  numbers  with  many  digits  more  readable,  an   additional 
feature  is  provided  in  the  format  specification  for  E  and  F  type 
conversion.   Either  the  F  or  the  E  specification  as  described 
previously  may  be  preceded  by  mK. 

The  effect  of  mK  is  to  cause  blocks  of  m  digits  to  the  right 
of  the  decimal  point  to  be  separated  by  one  blank.   The  separation 
will  occur  on  both  E  and  F  type,  and  it  will  continue  until  reset 
by  another  mK  appearing  in  some  specification. 

When  this  feature  is  used,  the  column  width  w  must  be  increased 
enough  to  allow  j  extra  spaces . 

Let  —  =  k  +  R  ,  k  an  integer  and  o  <L   R  <  m. 

m  - 

o  if  m=o. 
Then  j  =  -{    k  if  R>o. 

k-1  if  R=o. 


Example ;   Given  Xo,DX,N,K 

arc  tan  (Xi)  ,  >: 

i  =  1, 2, . . . ,N 


Xi 
1.   ComputelTe    and  arc  tan  (Xi) ,  Xi  =  Xo+(i-l)DX 


Xi 
2„   Print  a  table  of  Xi  ,  TTe    ,  arc  tan  (Xi) 

3.   If  K  ^   ±o  then  use  MODIFY  to  reset  the  package  so  that 
|k|  bits  are  used  for  the  power  of  two  exponent. 

4o   If  K<   -o  give  an  octal  dump  of  arc  tan  (Xi)  ,  i  =  1,2,...,N. 

5.  If  Xo  =  DX  =  o  then  terminate  the  job  by  CALL  EXIT. 

6.  MPP  to  be  used  is  a  relocatable  subprogram  which  operates 
under  FORTRAN  II  under  IBSYS. 
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$EXECUTE        FORTRAN 

$ID    THECODER*001/63/003$EXAMPLE  OF  MPP  USE 

*  XEQ 

*  FAP 

* EXAMPLE  OF  USE  OF  MPP 


COUNT 

90 

EXTERN 

ENTRY, C LA ,S TO 

, DPFBCD , DPEXP, MPY, DPATN 

EXTERN 

DPLOAD, MODIFY 

,  TEST 

, EXIT , ADD , ETRA , DPEXI 

EXTERN 

D(BDC) ,D(FIL) 

rDPBDC3 

START 

TSX 

DPLOAD, 4 

READ  N  AND  K 

TXL 

N,0,K 

DONT  STOP  ON  ERROR 

CLA 

K 

TZE 

NOMOD 

MOD 

STA 

*+2 

TSX 

MODIFY, 4 

MODIFY  MPP  TO 

PZE 

** 

USE  N  BITS. 

NOMOD 

TSX 
TXL 

DPLOAD, 4 
XO,0,DX+2 

READ  XO  AND  DX 

TSX 

TEST, 4 

CHECK  CARD  ERRORS 

TMI 

EXIT 

YES EXIT 

CLA 

DX 

NO 

TNZ 

AA 

CLA 

XO 

TNZ 

AA 

CALL 

EXIT 

DX=AO=0 — SO  EXIT 

AA 

LXA 
AXT 

N,l 
0,2 

TSX 

ENTRY , 4 
DPFBCD,0,PII 
STO,0,PI 
DPEXIT 

ENTER  INTERPRETIVE 
REPLACE  BCD  PI 
BY  M.P.N.  PI 
LEAVE  INTERPRETIVE 

BB 

TSX 

TXI 
TIX 

ENTRY, 4 

CLA,0,XO 

ST0,2,ARGT 

DPEXP 

MPY,0,PI 

ST0,2,EXPT 

CLA,0,XO 

DPATN 

ST0,2,ATNT 

CLA,0,XO 

ADD , 0 , DX 

STO,0,XO 

ETRA,0,*+1 

*+l,2,-3 

BB , 1 , 1 

ENTER  X(I)  IN  TABLE 
EXP(X(I)) 

ENTER  EXP  IN  TABLE 

ATAN(X(I) ) 

ENTER  ATAN  IN  TABLE 

INCREMENT 

LEAVE  INTERPRETIVE 

• 

PRINT  THE  TABLE 

TSX 

HEAD, 4 

START  HEAD  PRINT 

TRA 

D(FIL) 

END  PRINT 
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TSX 

TRA 
HEAD   TRA 

BCI 
TFORM  TRA 

BCI 
DHEAD  TRA 

BCI 
DFORM  TRA 

BCI 


TFORM, 4 

N,l 

0,2 

CLA,2,ARGT 

CLA,2,EXPT 

CIA  , 2 , ATNT 

*+l,2,-3 

LOOP, 1,1 

D(FIL) 


START  TABLE  PRINT 


PICK  UP  X 

PICK  UP  PI* EXP (X) 

PICK  UP  ATAN(X) 


GO  BACK 
END  TABLE 


PRINT 


CHECK  FOR  OCTAL  DUMP 

K 

START 

0,2 

0,2 

*+l,2,ATNT-l 

*+4,2 

DHEAD , 4 

10,0,10 

DFORM, 4 

ATNT,0,** 

START 

D(BDC) 


NO  DUMP  IF  PLUS 


HEAD  THE  DUMP 


PRINT  THE  DUMP 

**=ATNT+N-1 

GO  BACK  FOR  NEW  DATA 

PRINTS  HEADING  ON  NEW  PAGE 

8 , 1C1, 12X1CX, 15X, 9CPI*EXP (X) , 20X, 10CARC  TAN (X) // 

D(BDC) 

6,1C  3)1F2.15,5K1S1P5(1E17.31,F17.27 

DPBDC3 

9,1C011X3CL0C5X,4C(. . .12X10CARC  TAN(X) ,12X,4C. . .) 

DPBDC3 

5,3S8)2P1012.32,2012.15 

6, D3. 14159265358979323846264338328 

1 

1 

3 

3 

3 

3*2000 

3*2000 

3*2000 

(MPP  PACKAGE  GOES  HERE) 


D-10 


5,0 

8,-11 

0,0 


N  AND  K 
XO  AND  DX 
N  AND  K 
XO  AND  DX 
END  JOB 
END  JOB 
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The  above  program  and  data  would  cause  two  pages  of  output 
as  follows. 


First  page 
X 


PI*EXP(X) 


ARC  TAN(X) 


■10.00  (-4)  1.42628  08581  53150  16 

-5.00  (-2)  2.11678  84792  60429  67 

0.00  (0)  3.14159  26535  89793  24 

5.00  (2)  4.66253  69033  27078  08 

10.00  (4)  6.91981  83125  51164  68 


■1.47112    76743    03734    59 

■1.37340    07669    45015    86 

0.00000  00000  00000  00 

1.37340  07669  45015  86 

1.47112  76743  03734  59 


Second  paqe 

X 

PI*EXP(X) 

ARC  TAN 

(X) 

-90.00 

(-39) 

2.57422 

49862  95062 

81 

-1.55968  56728 

97289 

15 

-50.00 

(-22) 

6.05934 

63529  75874 

74 

-1.55079  89928 

21746 

09 

-10.00 

(-4) 

1.42628 

08581  53150 

16 

-1.47112  76743 

03734 

59 

30.00 

(13) 

3.35725 

50038  09131 

03 

1.53747  53309 

16649 

42 

70.00 

(30) 

7.90248 

36491  15328 

54 

1.55651  15842 

07499 

99 

110.00 

(48) 

1.86012 

82224  22199 

02 

1.56170  56681 

29836 

80 

150.00 

(65) 

4.37846 

77798  59209 

79 

1.56412  97588 

91028 

39 

190.00 

(83) 

1.03062 

68066  99454 

58 

1.56553  32174 

97301 

24 

LOC 

(... 

ARC  TAN(X) 

...) 

(51474) 

707507436664 

610642266434 

433612752001 

(51477) 

706401123263 

533216054600 

722566142001 

(51502) 

674233645405 

600051602725 

576063242001 

(51505) 

304627767342 

306340135451 

171727346001 

(51510) 

307167426070 

020602533113 

015270442001 

(51513) 

307713742512 

141612203176 

027716456001 

(51516) 

310152635504 

214501465665 

261677306001 

(51521) 

310306621707 

317331605432 

265347206001 
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CLA1 

CLS 

GET3 

LDQ 

SSM 

SSP 

SUB 

Reference  List  of  Entries  In  MPP 

Non-control  Pseudo-operations 

ADD  CACMQ  CHS  CLA 

CLS1  DVH  DVH2N  GET1 

MPY  MPY2N  PUT1  PUT3 

STO  STOl  STQ  STZ 

Control  Pseudo-operations 

DPEXIT   ETMI     ETRA     ETZE     TMI      TNZ 
TPL      TXI      TZE 

Elementary  Function  Pseudo-operations 

ATOB     DPACOS   DPASIN   DPATN    DPCOS    DPCOSH 
DPEXP    DPLNX    DPS IN    DPSINH   DPSQRT 

Data  Transformation  Pseudo-operations 

DPFBCD   DPFLTF   DPFLTI   DPTOM 

Miscellaneous  Pseudo-operations 

LXA      LXD      SXA      SXD      TIX      TNX 
TXH      TXL 

Control  Subroutines 

ADJUST   ENTRY    MODIFY   NOSCAN   SETERS   TEST 

Input  Subroutines 

DPLOAD   RTX 

Output  Subroutines 

DPBDC3   D(BDC)   D (FIL)   WOT 

Mathematical  Subroutines 

EQS 
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